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We evaluate the tree level four fermion string amplitudes in the TeV string mass scale models 
with intersecting Z)6-branes. The coefficient functions of contact interactions subsuming the con- 
tributions of string Regge resonance and winding mode excitations are obtained by subtracting out 
the contributions from the string massless and massive momentum modes. Numerical applications 
are developed for the Standard Model like solution of Cremades, Ibanez, and Marchesano for a 
toroidal orientifold with four intersecting D6-brane stacks. The chirality conserving contact inter- 
actions of the quarks and leptons are considered in applications to high energy collider and flavor 
changing neutral current phenomenology. The two main free parameters consist of the string and 
compactification mass scales, rus and Ale- Useful constraints on these parameters are derived from 
predictions for the Bhabha scattering differential cross section and for the observables associated 
to the mass shifts of the neutral meson systems K — K, B — B, D — D and the lepton number 
violating three-body leptonic decays of the charged leptons /i and t. 
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I. INTRODUCTION 

The consideration of Dirichlet branes has led in recent years to remarkable advances in the particle physics model 
building. This has made possible the construction of wide classes of Standard Model realizations for type // superstring 
theories using branes which extend along the flat spatial dimensions of M4 and wrap around cycles of the internal 
space manifold. The existing two approaches employing configurations of multiple type lib branes located near 
orbifold singularities [1-3] and type I la branes intersecting at angles [4], which we designate henceforth for lack of 
better names as setups of branes within branes and intersecting branes, respectively, are well-documented by now, 
thanks to the reviews in [5] and [6-8] . The two most characteristic features of these string constructions reside in the 
wider freedom in choosing the string theory mass scale and in the occurrence of localized chiral fermions in the string 
spectrum. 

Because they are amenable to experimental tests based on data from high-energy colliders and limits for rare 
processes, the TeV string mass scale models [9, 10] are clearly those with the greatest impact on phenomenology. In 
developing the string theory machinery for the physics from extra dimensions, one is especially encouraged by the 
applications using orbifold field theories in higher dimensional spacetimes with matter fermions which move in the 
bulk or are localized inside thin [11, 12] or thick [13-16] domain wall branes. The theoretical and experimental aspects 
of collider studies of physics from extra dimensions are reviewed in [17-20] and [21]. The string theory framework 
embodies a higher degree of consistency in comparison to the field theory framework, and is also more economical 
thanks to a small parameter space restricted to the fundamental tension and coupling constant string parameters, 
ml — 1/a' and gs, along with adjustable parameters associated with the compactification and infrared cutoff mass 
scales. Several examples that qualify as TeV scale string models have been constructed within the branes within 
branes [22] and intersecting branes [23] approaches. 

One important motivation for the interest in TeV scale string models is to gain insight on the hierarchy between the 
contributions from the various string excitations. In particular, one wishes to understand how the exchange of string 
Regge and winding modes compares with that of the gravitational Kaluza-Klein (KK) modes. This has a practical 
importance since, unlike the contributions the latter ones can be well described, in principle, within the familiar field 
theory framework. 

The early string inspired studies were focused on models using single brane configurations [24, 25]. These have the 
characteristic property that the open string Regge resonance modes contribute at the tree (disk surface) level while the 
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closed string modes contribute at the one-loop (cylinder surface) level of the string perturbation theory on the world 
sheet. (That the exchange of closed string modes can be viewed as a loop effect follows from the world sheet duality 
linking the ultraviolet and infrared regimes.) The strong constraints from the J\f = 4 supersymmetry preserved by the 
extremal branes of type II supergravity which saturate the BPS (Bogomolny-Prasad-Sommerfeld) bound, restrict the 
low-energy tree level contributions in single brane models to local operators of dimension V > 8. Using the orbifolding 
mechanism, say, by placing the Z?-brane at an orbifold fixed point in order to reduce the number of preserved 4-d 
supercharges and produce a chiral open string spectrum, is indeed helpful towards obtaining semirealistic models but 
does not significantly modify the structure of string amplitudes. This is often referred as the inheritance property of 
orbifolds. Another characteristic of the single Z?-brane models is the insensitivity of predictions with respect to the 
structure of the internal space manifold. Here again the orbifold constraints on the Chan-Paton (CP) gauge factors 
greatly improve the description by ensuring that massless poles in channels with exotic quantum numbers cancel out. 
Useful applications to neutrino-nuclcon elastic scattering at ultra-high energies [26, 27] and to the two-body reactions 
at high-energy colliders [28-30] have been developed along these lines in terms of single D-brane models where mg 
and the CP factors are treated as free parameters. Also, building up on the initial studies of the single photon-|-jet 
signal in the high energy hadronic collider reaction, p + p ^ 1 + j, recent works discuss the dijet signal [31] in the 
reaction, p + p ^ j +j, along with signals from the various processes at the LHC model [32] which could usefully test 
the single D-brane models of TeV string scale.. 

The multiple brane setups bring a crucial novel feature in this discussion through the presence of localized massless 
fermions whose contact interactions are not restricted to dimensions D > 8. This point was first recognized by 
Antoniadis et al., [33] in the context of TeV scale string models with Dp/ D{p+4)-hra.ncs. Finite tree level contributions 
were indeed found for the dimension V = 6 local operators coupling four fermions of which at least a single pair belongs 
to the non-diagonal open string sectors {p,p+4:) + {p+4:,p). That the open string fermions localized at the intersection 
of _D-branes behave in much the same way as the twisted modes of closed strings [34, 35] was realized by several authors 
in the context of multiple brane [36 38] and intersecting brane [39 44] models. The operator algebra approach for 
the superconformal field theory on the world sheet can be used to calculate the tree level string amplitudes. We also 
note as a side remark that the tools developed in recent years for the calculation of scattering amplitudes in gauge 
theories [45] and string theories [46] should encourage pursuing applications for general n-point amplitudes at higher 
loop orders. 

So far, the collider studies of TeV string scale models have been mostly focused on setups with single branes [24, 25] 
and branes within branes [33, 47, 48], as already said above. No comparable applications exist for the intersecting 
branes models. Regarding the flavor physics, however, several studies have been devoted to both the branes within 
branes [22] and the intersecting branes [49, 50] models. In parallel, a wide interest was aroused by the field theory 
models in flat extra dimensions with thin branes [51, 52] and thick branes [53-57] or in warped spacetimes [58, 59]. 
The reader should be warned that the quoted references represent a tiny fraction of the literature on this subject. 

In the present work we wish to pursue the discussion of the tree level contact interactions for four localized fermions 
in intersecting brane models with the view to confront the predictions against experimental data for colliders and 
flavor changing neutral current processes. Although the assumption of a low string mass scale is naturally paired 
with that of large extra space dimensions, the requirement that the string theory remains weakly coupled turns out 
to restrict the ratio of string to compactification scales. m,sr = nis/Mc, to a relatively narrow interval of 0(1). This 
circumstance has motivated us in taking the contributions from world sheet instantons into account. We evaluate 
the tree level string amplitudes by integrating the vacuum world sheet correlators over the moduli space of the disk 
surface with two pairs of massless fermion vertex operators inserted on the boundary. Since the tools for calculating 
open string amplitudes in intersecting brane models [36-39, 41, 43] are well documented by now, we shall present 
very briefly the main formulas before proceeding to our main goal. The discussion will rely heavily on our previous 
work [44]. 

We develop concrete calculations for the Standard Model solution obtained by Cremades et al., [60, 61] for a toroidal 
orientifold with four D6-branes, using specifically the related family of solutions presented by Kokorelis [62]. This is 
conceived as a local model (or premodel) described by a classical configuration of intersecting D6-branes decoupled 
from the gravitational interactions and the geometric moduli fields. It is encouraging that other searches of solutions 
for oricntifolds with intersecting branes also select families of small size, as illustrated by the study focused on the 
Z2 X Z2 orbifold models [63] realizing the supersymmetric Pati-Salam model with a hidden sector, and the statistical 
studies of the open string landscape of vacua for the Z2 x Z2 [64, 65] and Zq-jj [66] orbifolds realizing the minimal 
supersymmetric standard model and the Pati-Salam model with hidden sectors. 

The proper identification of physics from extra dimensions presupposes that one can combine the new physics 
contributions with those from the Standard Model interactions in a consistent way. This condition is especially 
critical for the string theory applications where a satisfactory implementation of the elcctroweak symmetry breaking 
is not yet available. Rather than pursuing a fuU-fiedged calculation, we shall adopt here a phenomenologically minded 
approach, similar to that used in [33] . This consists in separating out by hand in the low-energy expansion of string 



3 



amplitudes the contributions from the string massless and momentum modes so as to access the contact interactions 
which subsume the contributions from the string Regge and winding excitations. 

The exchange of massive modes from extra dimensions can also induce flavor changing interactions among fermions 
of different flavors which sit at points finite distances apart along the extra dimentions. These effects come on top of 
the flavor mixing effects generated during the electroweak gauge symmetry breaking by the trilinear Yukawa couplings 
of fermions to Higgs bosons. We focus here on a restricted set of hadronic and leptonic flavor observables believed to 
be among the most sensitive ones. To simplify calculations, we introduce certain assumptions on the flavor structure 
of the four fermion amplitudes which lead to an approximate factorization of the direct and indirect flavor changing 
effects. 

The outline of the present work is as follows. Building up on our previous work [44], we present in Section II 
the tree level four fermion string amplitudes for the high energy processes of fermion-antifermion annihilation into 
fermion- antifermion pairs and fermion pair scattering, / + /—>/' + /' and f + f f + f ■ We next consider 
an approximate construction of the contact interactions between pairs of quarks and/or leptons produced by the 
decoupling of string excitations. Finally, specializing to the Standard Model solution of Cremades et al., [60, 61], 
we present numerical results for the chirality conserving contact interactions as a function of the string and com- 
pactiflcation mass scales and the parameters describing the separation of intersection points. The corrections to the 
Standard Model contributions are studied over the admissible parameter space for the string and compactification 
mass scales. For comprehensiveness, we provide in Appendix A a brief review of the intersecting D6-brane models 
putting a special emphasis on the topics relating to the parameterization of the branes intersection points and the 
Chan-Paton gauge factors which have been lightly addressed so far. The discussion of tree level string amplitudes is 
complemented in Appendix B by a review encompassing both the intersecting £)6-brane and DS/DT-brane models 
aimed at the two-body processes of fermion-antifermion annihilation into pairs of gauge bosons and of gauge boson 
scattering, / + / ^ 7 + 7 and 7 + 7-^7 + 7. 

In Section III, we discuss the implications from the indirect high energy collider tests with a special focus on the 
Bhabha scattering differential cross section. In Section IV, we examine the contributions from the flavor dependent four 
fermion contact interactions to the hadronic and leptonic flavor changing observables associated to the mass splitting 
of quark-antiquark neutral mesons and the lepton number violating three-body leptonic decays of the charged leptons. 
For all the above applications, we compare our predictions with experimental data in order to infer lower bounds on 
the string mass scale at a fixed ratio of the string to compactification mass scales. 

II. TREE LEVEL STRING AMPLITUDES IN INTERSECTING BRANE MODELS 

We calculate the tree level open string amplitudes for four massless fermion modes localized at the intersection of 
£'6-branes in toroidal orientifold models. After quoting in Subsec. II A the general formula for the string amplitudes 
and discussing its low-energy representation as infinite sums of pole terms and the subtraction prescription proposed 
to construct the contact interactions, we specialize in Subsec. II B to the Standard Model vacuum solution of Cremades 
et al., [60, 61] and present in Subsec. II C numerical predictions for the contact interactions. Several notations are 
specified in Appendix A which provides a brief review of intersecting D6-branes models. All calculations are performed 
with the space-time metric signature, ( — using units where, a' = l/m^ = 1, except on certain occasions where 
a' will be reinstated. 



A. Four fermion string amplitudes and contact interactions 

The open string amplitudes for four fermions localized at points of the internal manifold can be calculated most 
conveniently by means of the superconformal field theory on the world sheet. The basic tools were initially developed 

for the closed string orbifolds [34, 35] and refined in several subsequent works (see [68] , for instance) . The application 
to the 'twisted' or non-diagonal modes of open string sectors was discussed later for the case of branes within 
branes [33, 36-38] and of intersecting branes [39, 41, 44, 49]. 

1. String amplitudes of localized fermions 

The general configuration of quantum numbers for the four fermion processes consist of two incoming conjugate 
fermion pairs, fi{k\) + ^(fe) + /^(fcs) + fi{ki), localized at the four intersection points, Xi, Xj, Xf., Xi e T^, of 

the four D6-brane pairs, {D,A), (A, B) and {B,C), {C,D), intersecting at the angles, and ^0 ^ . The tree level 
open string amplitudes are obtained from the correlators of vertex operators inserted at points xi, X2, X3, X4 on the 
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disk surface boundary by integrating over the disk moduli space of the punctured disk. Using the invariance under 
the Mobius group to set, Xi =0, X2 = x, 3:3 = 1, X4 = X ^ 00, one can write the resulting formula as [44] 

A'f4. = A{V_0^(^D,A),iMi^iWe,iA,B)JMi^-2W-e',(B,C),kMi^3)V0'^^^C,D),lMi^'i))/['>^^^ 



= C |^<Si234^1234Vl234(s, t) + <Si324'7"i324Vl324(w, t) + <Si243^1243Vl243(s, u) j , (II.l) 

where we use the notations: 

Jo J 

"71234 = 7l234 +T4321, T1234 = 7'r-(AiA2A3A4), <Si234 = -«Sl432 = (uf 7°7/:tW2)(M3 7°7''W4), 

S = -(fci + k2f, t = -(fc2 + ksf, U = -(fci + ksf. (11.2) 

The string amplitude A!f4, in Eq. (II.l) is built from three reduced (partial) amplitudes associated to the cyclically 

inequivalent permutations of the insertion points. The second and third terms are obtained from the first term 
associated with the reference configuration 1234 by substituting the mode labels 2^3 and 3^4, and modifying 
the interval of the cc-integral from x G [0, 1] to cc S [l,oo] and x S [— oo,0]. Each partial amplitude decomposes into 
a pair of amplitudes associated to the direct and reverse orientation permutations of the labels, corresponding to the 
substitutions, 2 4 and x ^ (1 — x). The requirement that the world sheet boundary is embedded in the Tj on closed 
four-polygons with sides along the brancs ABC'D. can be satisfied, in general, only by a single partial amplitiide, the 
conflict with the target space embedding forcing the other two to vanish. In our present notational conventions, only 
the reference term, 71234V1234, survives, while the other two terms, 71324V1324, ^243Vi243, cancel out. 

The factorization of the residues of the massless pole terms from exchange of gauge bosons between fermion pairs 
into products of three point current vertices determines the normalization factor as, C = 2Trgsa' . Using the familiar 
results [69] for the D-branes of type II string theories yields the formula in Eq. (II.l) expressing C in terms of the 
gauge coupling constant (/^ of the 4-d gauge theory on the _D6;j-brane and the volume of the three-cycle that it 
wraps. The same formula applies to each factor of the complete gauge group. The orientifold symmetry is taken into 
account by the factor K/j, which is assigned the value K/^ = 1 or Kfj_ = 2 when the brane /i is distinct or coincides 
with its mirror image, corresponding to the cases with U{N) and extended SO{2N) or USp{2N) gauge symmetries, 
respectively. It is important to realize that the relation between the string theory gauge coupling constants, 5^ , and 
their field theory counterparts which we denote momentarily by g^^*, also depends on the way in which the analogous 
field theory model is constructed. An illustrative discussion of the model dependence is presented in [43]. For the 
moment we express this relationship by the proportionality relation, rjgjj* = 3^, involving the real parameter 77. 

The factors Ii{x) and Z^^ = Xlc/ ^"^^^ each in Eq. (II.l) designate the quantum (oscillator) and zero mode 
world sheet instanton contributions to the correlator of coordinate twist fields in the complex plane of Tf. With 
the choice of independent pair of cycles, Ca — {xitX2), Cb = (2^2, 2:3), surrounding the insertion points along the 
world sheet boundary, xi, ■ ■ ■ , x^, which map to the intersection points labeled i,j,k,l in Tf , the summations in 
the classical partition function run over the large lattice generated by the one-cycles and wrapped by 
the 2?6A/-D6s-branes in T/ . Going through a full circle around the cycles Ca, Cb, induces the coordinate fields 
monodromies 

V2/^c^X = 27r(l - e^''''^)vA, ^Ac^X = 27r(l - e^"^)^;^, (II.3) 

where 

VA = PaLa + Sfj = {pa + efj)LA + , vb = PbLb + Sf^ = {pb + t%)LB + df^, (11-4) 

with Pa, Pb & Z denoting the winding numbers. The 2-d large lattices generated by the brane pairs in the complex 
planes of Tf are displaced from the origin by the shifts separating the branes intersection points, = e^LA+d^, 6^ = 
€^ Ls+d^ , [^{L\ gd^'^) = 0] where the real parameters e"^, and rf^, d,^ stand for the longitudinal and transverse 
components of the shift vectors relative to branes A, B. Detailed formulas for the functions I{x) and can be 
found in our previous publication [44]. 

In the special case involving equal interbrane angles, 6^ = 6 ^ , which corresponds to the parallelogram DABC with 
D = B and A = C, the string amplitude simplifies to 

nl 



— C<Si234 



T12.34 / dxx-'-\l - a;)-*-i V 34(0;) 
•^0 cl 
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+T4321 / dxx-*-' (1 - X)-'-' V Zg^^.ix) 
•^0 cl 

where 

Zf^^^{x) ^ d('^'^\{TA)^C'^''^]{TB) = ^ g-'rsm7re7[|(pA+ef2)Li+df2| = -^:i;j^ + |(pB+ef3)Lij+df3|^T|^^]^ 

PA, PB&Z 

TA{x)=ism{'Ke')\L'^\^ ^^^~^^\ rB(ar)=isin(7r^^)|L^|2-^^(^, = 1 - ^; 1; a;), (II.6) 

with F{a,h;c;x) denoting the Hypergeometric function and ^?[^](t) the Jacobi Theta function. We have separated 
out the direct and reverse permutation terms in the quartic order trace factor, T\2Zi- Verifying the equahty of the 
factors multiplying Ti2;m and T4321 provides a useful check on calculations. 

The string amplitudes for four fermions localized at the intersections of Dp/D{p + 4)-branes are derived [33] by a 
similar method to that used for intersecting branes. The resulting formulas are detailed in Appendix B along with 
the similar results for the two-body processes, / + / — > 7 + 7 and 7 + 7 — > 7 + 7. 



sin nO 
F(x)F(l-x) 



(11.5) 



2. Low-energy representations 

The low-energy limit of compactified string theories is described by means of series expansions in powers of y/s / rris 
and (m^r)"^ = where y/s stands for the energy variable and r = 1/Mc for the characteristic compactification 
radius parameter. In the flat space limit, r ^ 00, both the partition function and fermion localization factors in the 
.r-integral of Eq. (II. 1) can be set to unity, Z^i = e"'^':'^^^ 1, Ii{x) — > 1. Ignoring momentarily the gauge 
factor, one can express the string amplitude in this limit by the formula 



B{-s,-t) = —S{s,t) = / dxx-'-'^il 
St Jo 



x)-'-\ 



r(i-s)r(i-t) r(-.)r(-t) _^ >t-i)---{-t-n) 

where wo have exhibited in the second line the representation of Eulcr Beta function in terms of an infinite series of 
polos from s-channel exchange of the string Roggo resonance of masses, — tittz^, [7^ Z^. At finite r, applying 
the familiar method of analytic continuation past poles to the x-integral [69] with the factors and I{x) included, 
produces infinite series of s-channel and t-channel pole terms at the squared masses of the open string momentum and 
winding modes. Since for the equal angle case associated with the brane configuration DABC with C = A, B = D, 
the argument variable ta{x) — *■ at small x, one must carry beforehand the modular transformation on the Theta 
function with modular argument ta{x). At small 1 — x, the same applies to the Theta function with modular argument 
tb{x). 

selecting the regions near x = in the representation of Eq. (II. 6) yields the low-energy expansion of the string 
amplitude in the open string compactification modes 

A> r^r<Z / -'1234 ih ^AB e PB 23 j^3^^ 1 1/ OgC ^ \ Ctt »^ 



where 



MZ^^b,, = sin^(7r^/)(PA + 4/f\L'A? + In^AB = 2^(1) - V(^') - ^(1 - O'), (II.9) 

and '4'{z) is the PolyGamma function. The squared masses, M\ consist of string momentum modes of the sector 
{B,B) along the compact directions wrapped by the brane B, and string winding modes of the sector (^4,^1) for 
strings stretched transversally to the cycle Lb- Expressing the normalization factor C = 2Trgs in the s- and t- 
channel polo terms in terms of the corresponding parameters of branes B and C, respectively, then comparison 

with the amplitude of the analogous field theory leads to the identifications, ^^Ti23a = (.9^ )^ Sa(^")i2(^")34 ^.nd 
2^ = (Sc )^ So(-^")23(^'')i4' where a labels the Lie algebra generators of the gauge groups on branes B, C. The 
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extra factor 2 in the denominator accounts for the field theory normahzation convention used for the field theory 
gauge coupling constant, corresponding to the gauge current vertex, < J^(0) >= gAT°'^x, [Trace{TaTi,) = ■^Sab]- 

The residues Sj^g ^''^ of the poles at s = M\ g (and the similar residues of the poles at t = represent the 

squares of the form factors for the three point couplings of fermion pairs localized at the D6a/DQb brane intersections 
with the open string states of mass Ma,b from the sectors {A, A) and (B,B). As the interbrane angle runs over the 
defining interval, 6^ G [0,1] in the Tf tori, ln6^g asymptotes to cot^nO^) +oo at the interval end points while 
reaching the minimum value, InS^^^g ~ 2.77, at the mid-point 0^ = ^. The divergence of ln(5j^^ at 6^^ = and 1 
implies then the absence of contributions from the exchange of string compactification modes for parallel branes, as 

— /2 

expected by virtue of the momentum conservation. The form factor, J^ab{pb) = ^ab ^'^ ' representing the cost for 
a fermion particle to absorb the gauge boson momentum ps, arises as a consequence of the D- brane fuzziness caused 
by the finite spatial extension of string modes. Its configuration space representation, J^abIv), can be evaluated by 
writing the three point coupling of fermion pairs to the momentum modes of the DGs-brane gauge connection field, 
A^{x,y), as a Fourier integral over the cycle of radius Lb 



£ " dyA^{x,y)J^AB{y-yn^^^^J2^PB[A^^'"'^{x)£ rfy cos(^)^^B(y - 2/f ) + • • •] 

V ' ^' PB 

with the resulting formula 

TAB{y) ^ \l , e"^i'^Bi-.n.,3^ [yS = 27ref (11.11) 

V nmOAB 

where y parameterizes the points on the Lb cycle which has been described here by the orbifold jZi of length 
7r|I/B|; yf denotes the position of the fermion mode; iVp^ = 1 for 7^ and N^^ = 1/a/2 for — 0; and the 
central dots stand for the sine stationary modes. The above approximate formula for TAB^y) becomes exact in the 
large radius limit, \Lb\ — > 00. The exchange of massive KK gauge bosons contributes to the four fermion amplitude 
a sum of pole terms with residue factors, ^'^'^VB^ju^ ^^I\^b\ ^ ^j^^j.^ ^b^ denotes the relative distance along the brane 
B, as expected by comparison with Eq. (II. 9). A similar analysis holds for the momentum modes associated with 
the brane C . By contrast, there is no field theory interpretation for the form factors accompanying the coupling of 
localized fermions to the open string winding modes. 

To illustrate further how the form factor originates within the field theory framework, we consider the toy-like 5-d 
f/(l) gauge theory with the fifth dimension compactified along the orbifold segment, y e [0, 7r|I/B|], assuming that the 
chiral fermions are trapped near the boundaries by some soliton kink solution involving a scalar field coupled to the 
fermions. Making use of the Gaussian ansatz for the normalizable zero mode wave function of a fermion localized at 

i,f\y) = 7Ve-(^-^?)V(2-^-^), [ dv^T\y)^f\y) = 1, AT (4-)'] (11.12) 

Jo 772(7 

where the normalization integral determining N has been evaluated in the limit Lb —> 00, we infer the gauge vertex 
coupling 

PB ^ ' ^' PB I S| 

Comparison with the form factor J^abIpb) in Eq. (11.11) allows us to identify the half-width parameter as, a ~ 
V21ni5yiB/7r. The above result for the form factor of fermion modes localized at intersecting branes agrees with that 
derived [33] for the fermion modes of the non-diagonal sectors of Z)3/i)7-branes, where 6^ = \. (The half- width 

parameter in [33], which we distinguish here by the suffix label ABL, is related to ours as, (jabl = cr7r/\/2. For 
later reference, we note that our half- width parameter a identifies with the parameter denoted a in [56].) It is of 
interest to note that a similar form factor also arises in the string amplitude for emission of a massive graviton mode 
G{Mci) in the two-body reaction, f + f ^7-1- G{Mci), with the characteristic dependence on the mass Md of the 
closed string mode [25], Tg{Mci) = e~ ''^*^^)^<!'/'"= . Both the present string form factor and the one derived above, 
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Top{Mop) = e"*^"? should be distinguished from the form factor arising from the quantum fluctuations 
of branes [70-72]. Using a quantum field theory treatment of the Nambu-Goto action for domain walls with the 
transverse coordinate y promoted to the would-be light Nambu-Goldstone scalar field associated to the spontaneous 
breaking of translational invariance, Bando et al., [70] found that the coupling to momentum modes in the effective 
4-d action is accompanied by the recoil form factor 

J-rec{M)-e , [A ^^^^^^^ 4^4 ^^^^ j^^pi l^^-^^J 

where Mp denotes the ultraviolet mass cut-off for the field theory on the domain wall whose tension parameter has 
been identified above to that of the Dp-brane, rop = [m1/{2-KgDp)\'^ , with gauge coupling constant gop- The specific 
case of a soliton domain wall is discussed in [71] . From the comparison of the ratio of half- width parameters for the 

form factors associated to the brane fuzziness and recoil, : cr^ = ^^f : 4^ — -rr^ '■ 1, where we have assumed 

■ 4m^ lino ' 

in the second stage, Mp — rUg, we conclude that the suppression effect from brane recoil is parametrically weaker 
than that from the string finite size. This property was previously discussed in [25]. 



3. Contact interactions and flavor structure 

The four fermion string amplitude receives infrared divergent contributions from the a;-integral boundaries which 
correspond to the massless s, t channel pole terms from gauge bosons exchange. One way to regularize these harmless 
poles is to subtract out the small regions in the x-integral near the end points a; = 0, x = 1, while adding up the 
corresponding pole terms by hand, as described by the subtraction, A—^ A — (— -I- ^|a;=o) ~ {—A^°^^ + ^|a;=i)- 
To account for the electroweak symmetry breaking, one can use the same prescription where the added pole terms 
correspond to the contributions from exchange of the physical gauge bosons with the observed finite values of the 
masses. For the four fermion coupling {fH^^fH)if'H'l^f'H') the electrically charge neutral channel for the 7, Z 
gauge bosons, this is illustrated by the substitution, ^ ^ 7 + ^ ; where aH{f) denote the Z boson vertex 

couplings. A similar replacement holds for the i— channel poles. 

The contact interactions, subsuming the contributions from the massive string Regge and winding excitations, can 
be constructed in a similar way by subtracting also the pole terms associated to the string momentum excitations, as 
illustrated by the schematic formula 

•^'contact = -4' — •^'s,KK ~ •^t,KK- (11.15) 

The same subtraction procedure was previously used to define [33] the contact interactions in models with Dp/D{p+ 
4)-branes. For consistency, we remove the s— channel terms only for the configurations of intersection points with 
i = j ov = 0, and the t— channel terms only for those with fc = Z or = 0. No subtraction of momentum 
modes is needed in the cases, i ^ j ^ fc or ^ 0, ^ 0, where the thresholds for the momentum modes are 
separated by a finite gap corresponding to the mass terms, e^|i^|^. One motivation for excluding the KK towers 
from the contact interactions is simply that it is always possible to include separately their contributions through the 
familiar field theory treatment, suitably generalized by the inclusion of form factors. In the flat space limit, Lb — > oo, 
with the brane form factors set to unity, S~^^''^ — > 1, the contact interactions from KK modes are described by the 
approximate formula 

^^^^ - 2^/''''^'''' 1^ {py\LB\'-s) - ^ (n-2) '^i234Ti234, (11.16) 

where Sn = r(n/2) '^i^h n = D — A denoting the number of real dimensions of the wrapped cycle. 

In the kinematic regime where mass parameters and energies arc negligible in comparison to rus, the four fermion 
local couplings are dominated by the chirality conserving couplings of dimension V = 6. The corresponding terms in 
the Standard Model effective Lagrangian are represented in the left and right chirality basis of fermion field operators 
by the general structure 



Leff = Yll(^ff'{fL7Mfh^f'L) + G'ff'ihl^ml'Rl^fR) 

JJ' 
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where G^^f = tjhh' %/2 ? with H, H' = {L,R). We follow the notational conventions for the coefficients G^Jf 

commonly adopted in the context of composite models [73], with r]HH' denoting relative ± signs and ^{Jj^i the 
characteristic squared mass scale parameters. For later reference, we also quote the effective Lagrangian of use in 
leptonic collider applications 

^^^^= E Ti7S^r^Gf/'(e-H7^eH)(/H'7'^/H'), (n.l8) 

f H^H, A4h,{1 + Sef) 

where S^f = 1 if e = / and ^ = if e 7^ /. The left chirality basis for the fermion fields is related to the mixed left-right 
chirality basis by, L i?^^, thus implying the identity, (Li7''L2) = (-Ril'^-^i)) ^^'^ Fierz-Michel identities for the 
quartic order matrix elements, 

(Li7''L2)(^37/.i?4) = 2(Iii?4)(^3i2), {LlJ^L2){L37^.L4) = -{La''L4){Ls-f^L2), (11.19) 

where L, R denote c-number (commuting) Dirac spinors. The comparison with Eq. (II. 1) yields the explicit formula 
for the dimensionless coefficient functions {L^) defined by 

HH' _ mjK^ jjH' _^ rnlK^ _ \L^\ /-^ _i ._t_i-pr, 2sin7r0^ 1 ^ 

such that the string four fermion amplitude reads as, Afi = gfff < 5i234 >hh'< ^234 >ff' ■ For convenience, 
we shall employ in the sequel a similar notational convention for the string amplitudes for fermions of fixed chiralities. 

To make contact with the amplitudes of physical processes involving the mass eigenstate fields, we need to perform 
the familiar bi-unitary linear transformations linking the above gauge basis to the mass eigenvector basis, / — + 
Vpf, f Vr' r, [f = lA which read in the left-right chirality basis as, Jh ^ V^Vi/, [H = L,R\. The flavor 
mixing matrices, v[, V/j, are determined through the diagonalization of the fermion mass matrices in generation 
space, V^MfVp = (Mj)diag, but in a partial way since the Standard Model contributions from the quarks and 

leptons only depend on the products, Vckm = ^l^l ^^'i ^' ~ ^l^l 1 where the sufhx label CKM refers to the 
quarks Cabibbo-Kobayashi-Maskawa matrix. The effective Lagrangian of dimension "D = 6 in the vector spaces of the 
fermions generation and mass basis fields can now be expressed as 

Leff = J2I1 ^"m [ifiHrfjH){fU'l,f'iH')]fiav + H.c. 

fj' H.H- 

= E E [ifiHrfjHKfkH'lMUase + H.C, (11.21) 

where 

'^ij,kl = ,^2j^^9ij,kl - ^i'j',k'l'\^Hii'^Hjj')\^H'kk'^H'll')^ '^ijM = ^'^J^A [i-i--^^) 

by using the familiar tensorial notation for the flavor and mass bases coefScients in Eqs. (11.20) which are labeled 
by the same indices i,j,k,l G [1,2,3]. The 4-point couplings G^^i of localized modes axe subject to geometrical 
selection rules which are expressed in terms of the shift vectors, w^Aj associated to the embedding polygon with sides 
D, A, B, C in each Tf, by the conditions [67] 

w^A + Wab + w^BC + w^CA = ^ mod [w^BA e A^^/A^], (11.23) 
where we use notations deflned just below Eq. (A. 9). 

B. Standard Model realization with four branes 

We here specialize to the solution of Cremades et al., [60, 61] with the brane setup consisting of four Z)6-brane 

(baryon, left, right and lepton) stacks of size Na = 3, Nb = 1. Nc = 1, Nd = 1, supporting the extended Standard 
Model gauge symmetry group, U{3)a x SU{2)b x C/(l)c x U{l)d- The weak gauge group, USp{2)b ~ SU{2)b, identifies 
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with the enhanced gauge symmetry associated with the overlapping pair of mirror D6;,/D6;,'-branes, and the hyper- 
charge with the hnear combination of Abehan charges, Y = ^ — ^ — As noted by Kokorehs [62], the present 
brane setup belongs to the family of solutions described by the winding numbers and intersection angles listed in the 
following table whose members are labeled by the discrete parameters, p = (1, 1/3), e = ±1, e = ±1, [ee = 1] with 
/3j = 1 — 6i = [1, i], [i = 1, 2] in correspondence with the cases of orthogonal and flipped tori TI3. 



Brane (A^^) 




{nl,ml) {nl,ml) ejr''^'' 


Susy Charges 


Baryon (A^^ = 3) 


(1,0) 


(i,3pe/3i) (i,-3pe/32) (O,±0„,T^a) 


r(i),r(4) 


Left {Nb = l) 


(0,ee) 


(^,0) (0,-e) {eb,0,TOb) 


Vr(2),r-(4)/ 


Right (A-c = 1) 


(0,e) 


(0,-e) (^,0) {±e,,T0c,O) 


r-(3),r-(4) 


Lepton {Nd = 1) 


(1,0) 


(i,3pe/3i) (i,-3pe/32) {O,±0d,TOd) 


r(i),r(4) 



The massless spectrum of left chirality multiplets localized at the intersection points includes three generations of 
quarks and leptons. For the specific choice of brane angles characterized by a single vanishing angle and a pair 
of angles equal up to a sign, the various brane pairs preserve Af = 2 supersymmetry each, provided the complex 
structure parameters of TI3 satisfy the relation, /3iX2 = /32X3- We have indicated in the last two columns of the 
above table the brane-oricntifold angles and the spinor weights of the conserved supercharges, with the upper and 
lower signs corresponding to e = e = ±1. The finite intersection angles can then be expressed as, 0^ = 6'c = 5, Oa — 
Od = ^ tan~^(3p^e/3iX2)- Only when e = e = — 1 do all the four branes share the common spinor supercharge, 

''(4) = ( ), implying the existence of an unbroken Af = 1 supersymmetry in this case. The assignment of open 

string sectors and gauge group representations for the quarks and leptons and for the two Higgs bosons is displayed 
in the following table in correspondence with the gauge group SU{3)a x SU{2)i, x U{l)a x U{l)c x U{l)d- 



Mode 


q 








Hd 


Hu 


Brane 


(a, 6) + (a, 6') 


(c,a) (c',a) (d,6) + (d,6') 


{c',d) 


{d,c) 


{c,b) + {c,b') 


(c,6)t + (c,6')^ 


Irrep 


3(3, 2)1,0,0 


3(3,l)_i,i,o 3(3,l)_i,_i,o 3(1,2)0,0,1 


3(1,1)0,-1, 


_i 3(1,1)0,-1,1 


(1, 2)0,1,0 


(1, 2)0,-1,0 



The massless scalar modes for the pair of up and down Higg bosons arise as the hypermultiplet of the sector (c, b) = 

{c,b') with Af = 2 supersymmetry, due to the coincidence of the branes c and b along the first complex plane T^. 
The effective gauge field theory for this family of models is free from anomalies, despite the fact that the brane setup 
fails to satisfy the RR tadpole cancellation conditions. These can be satisfied, however, by including a hidden sector 
of distant brane stacks which do not intersect with the observable brane stacks. The quark generations [60, 61] are 
located in the three complex planes of Tf at the intersection points 



.(1) -k .(3) \ /pW(« 

1-36^^V' W^j(d^) 




where the three entries in the upper and lower arrays stand for the coordinates in the orthogonal reference 
frames, 3?(A-^)/(27rr() and 3(A^)/(27rr|), [/ = 1,2,3]. The intersection points for the lepton generations, 
-^^6^(0, ^c'^di^^)' -Pio^(z^'^), are described by similar formulas to those for the quark generations with Cg^' — > e;^\ e^^' — > 
ii^K The intersection points along the branes are labeled by the integer index k G [0, 1, 2] ~ [0, 1, —1], with the branes 
transverse distance described by the real parameters, e^^^j, e^^j, e^^j, e"^], rtc- (We deviate only in the defi- 

nition of the parameters in the third complex plane with the notations of [60, 61] which use the choice of unit of 
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length. 3(X-'^''^)/(67rr2^^). The relationship between us and them reads, = ''i{e]f')Lhem-) Numerical studies 

for subsets of these parameters have been reported [49, 50] in attemps to fit the quarks and leptons Yukawa coupling 
constant matrices. Directions to improve certain shortcomings of the predictions for the fermions mass matrices are 
reviewed in [74] . Two important features of the present family of models are that the intersection points are separated 
along the branes by distances of order, e ~ 1/3, and that all three intersection points for the electroweak doublet and 
singlet fermion modes are placed at points which lie at finite distances apart only in the single complex planes, 
and T|, respectively. 

The gauge matrices are constructed along the lines traced out in Appendix A. To the electroweak singlet leptons, 
ef = (c'.d), v'^ = {d,c), with multiplicities, Ic'd = —3, led = 3 and U{l)c x U{l)d charges, {Qc, Qd) = (—1)—!) and 



A 



(cd) 



with non-vanishing 



(+1,-1), we ascribe in the subspace of gauge quantum numbers (c, c?) the matrices Agc'^^ 
entries 7 and a, respectively, using the notations of Eq. (A. 10). Exphcitly, 7 = 1 for e'^ ~ (c',d)_i._i and a = 1 for 
v'^ (d, c)_i,i, with other entries set to zero. For the modes charged under the color group SU{'i)a, choosing the 
subspace (a, c), we ascribe to the electroweak singlet quarks, u'^'^ ~ (a, c), (f^^ ~ (a, c'), of charges (Qa, Qc) = (+1, Tl), 
the bifundamental representation matrices a|^"^' , A^"'^' with the non- vanishing entries for the column and array vectors, 
Oia, la, in the notations of Eq. (A. 10), transforming under the SU{i) group hmdamental representations. For the 
electroweak gauge group, U Sp{2) ~ SU (2), supported by the i?6b-brane with the U Sp projection, the two components 
of the electroweak doublet lepton mode, I = {v e)^ ~ (d, h), with charges, {Qd, Tf ) = (1, ±1), are ascribed the matrices, 
with the non- vanishing entries, 7 and a, using the notations of Eq. (A. 10). A similar construction holds for the 
two components of the electroweak doublet quark mode, q = (u, d), which are ascribed the matrices \q°'^\ with the 



non-vanishing entries, 
formulas 



and aa- The CP matrices of the quarks and leptons modes are then given by the explicit 



A„c — — ^ 




(11.25) 



where the normalization condition is of form. 



pair of branes vector subspaces associated to each mode. The quartic traces of CP matrices axe easily calculated from 
the above explicit representations. The sums of traces for the direct and reverse orientation terms, indicated below 
by the suffix label (rf + r) , are given by the formulas 

1 

2' 



and, for convenience, we have kept a record of the 



'^'^(d+r){X\e)\e)X\e)\e)) = 7'r(rf+r)(A|^,jA 



(e=)A(gc)A(ec)) 



1, ^^(d+r)(AL)A(ec)A(e)AL) = -, 



rrw+.)(Aj^^^A(gg)A|^)A(e)) = rr(d+r)(A(dc)A|^e)A(ec)Aj^e)) = Tr(^d+r){\uc)X\^c)X\^c)X{ec)) = ^5^0, 



^^(d+r)(A(q„)A|g^,A| )A(gc) 



[^■0/375 = ;r(^c«/3^75 + ^ab^-i^), q = u,d\. 



(11.26) 



Making use of the U{Nc) group identity, 2^^^ ^(T")q,^(T")^5 = {SasSfij — -^Sa/sSys), with the first and second 
terms inside parentheses being associated to the non-Abelian and Abelian group factors SU{Nc) and U{1) of U{Nc), 
one can cast the tensorial structure of the trace factor for quarks in the fundamental representation of SU {Nc) into 
the operator form 



"71234 = Tai3jS = -^{SafjSjS + SagSfj-y) = (T")q,/3(T'')^5 + 

^ 0=1 



+ 1 



2N^ 



(11.27) 
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where the traceless Lie algebra generators of SU{Nc) are normaUzcd as, Trace(TaTb) = ^5ab- (The right hand side 
is symmetric under the substitutions, f3 ^ S or a ^ "f.) The corresponding quartic trace for Icptons is 71234 = 1- 
Identifying the leading pole term in Eq. (II. 9) to the pole term in the analogous field theory, and using the low-energy 
limit with due account for the gauge factors in Eq. (11.26), one finds that the ratio of the string to field theory gauge 
coupling constants of the non-Abelian gauge factors should be set as, rj = \f2. 



C. Numerical results for contact interactions 



We start off by stating the main simplifications made in our numerical study. We choose to set all the radii 
parameters, r(, equal to a common radius parameter, denoted by r = 1/Mc. The relation between the string theory 
parameters, the wrapped three-cycles volume and gauge coupling constant parameters, gs, m^, La, Qa, \A. = a, b, c, rf] 
simplifies then to 

-«^=(5^)^' ['^>^=n^=n(-"+<v^)^ x'='f] (11.28) 

yAi'^M I I 

where Ka.c.d = 1, Ki, = 2. From the explicit formulas for the three-cycles volumes, we infer the relations between the 
branes gauge coupling constants 

Ca=Cd= [{p-^ + (3p/3iX2)')(p-' + (3p/?2X3)')]^ Cb = Cc = ^ 

Pi P2 

=^ 9a = gd, —5- = -^ = -T— • (11.29) 

9 c '-b P2X3 

In the case with A/" = 1 supersymmetry, /32X3 = /3iX2i one has, Cc = a-nd 1g^^ ~ g^^- Specializing momentarily 
to this supersymmetric case and using the proportionality relation between the string and field theory gauge coupling 
constants discussed after Eq. (II. 1), we can express the Standard Model gauge couphng constants 33, 52, 5i along 
with the electric charge and weak angle parameters, e and sin^^^, by the formulas valid at the string mass scale 

W3 = ga, Vg2 = gb, (Wl)"^ = ^ + ^9c^ + ^9d^ = ^ + 



2 3a~^ 1 

(.e)- = (..r)- + (.,2)- = + 2,-, sin^ 0^ ^ = ^-3/^ = 2(1 + ,g/(3.g)) -^"'^"^ 



Although the family of models under consideration has a parameter space of restricted size, to study the model 
dependence of predictions we found it convenient to introduce a reference set of natural values for the geometric 
parameters and consider small excursions in which the parameters are varied one by one. We define our reference set 
of parameters as, p = l,e = 1,? = = l,/32 = 1- 

The interbrane angle parameters in the three complex planes, 9j^^ = dl — [T^Sj^ = arctan(m^x^/^^)] a^^e grouped 
into two distinct sets with the entries in each set being equal up to permutations of the planes. For the reference 
set of parameters with = 1, we find 6^ (q) = O^e) = (0.50,0.60,0.89), and e\v') = 9^ [u") ~ 6^ {(F) = 9' {e") = 
(0.50,0.89,0.60). Let us also quote, in reference to the discussion near Eq.(II.ll), the numerical values assumed by 
the form factor parameter, In^^ = (2.77, 2.97, 10.1), for the above quoted values of the interbrane angles for 9^{q). 
Varying the tori complex structure parameters inside the range, G [5 '2], or changing from orthogonal to tilted 
tori, causes insignificant changes in the angles. For instance, the choice /3i = /?2 = | yields 9^{q) = (0.5,0.687,0.812). 

With the free parameters consisting of nis and m^r = m^/Mc, the string coupling constant is fixed in terms of 

the gauge coupling constants by, gs = —^^^^j^- For the string theory not to be driven to strong coupling in the 
decompactification limit at fixed rrig, the condition gs oc (mgr)^ < 1 restricts the radius parameter to m^r = 0(1). 
More precisely, setting tentatively in Wgr = ( ^g'i\cA\ ^ ^ ' ~ ^1 ^^"^ Standard Model gauge coupling constants to 

their observed values at the Z boson mass scale, gfimz) = 0.127, giifnz) = 0.425, gl{mz) = 1.44, we can express 
the conditions that the string theory is weakly coupled by the numerical results evaluated for the reference set of 
parameters, rusr ~ [0.95, 3.9]gl^^/r]'^^^. 

The weak angle depends sensitively on the geometric parameters. While setting ga — gb reproduces the favored 
value, sin^ 9w = |, using the reference set of parameters with x^ = I7 would yield instead, g1^/{g1/2) = Lb/ La — 
1/10 => sin^ ew ^ 3/46. However, as verified from the explicit formula, sin^ 9w = l/[2(l+2/3i(l+9/?Jx2)/(3xiX3))], 
one can always fit the weak angle by adjusting the geometric parameters, for instance, by setting, /3i = 1/2 with all 
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equal. The renormalization group analysis of the gauge coupling constants for the present model is known to be 
consistent with grand unification with the parameters adjusted at [75] rrig ~ Mx — 10^® GeV for m,g/Mc ~ (2 — 5). 
On the other hand, the qualitative study of the coupling constant unification for the class of minimal supersymmetric 
standard models selected by sampling over vacuum solutions of intersecting branes on orientifolds [76] indicates that 
the maximal allowed string scale may vary over a wide interval, provided one regards the branes gauge coupling 
constants as free parameters. In order to justify the TeV string scale scenario of interest to us, one must invoke string 
threshold corrections which produce an accelerated power law running between the compactification and string mass 
scales [77]. Although the small extent of the admissible interval for the momentum scale, Q G [Mc, nig], could render 
this option problematic, some freedom is still left in choosing the geometric parameters. To verify this statement we 
have pursued a qualitative study of the relation between the one-loop order running gauge coupling constants holding 
for the present model, 

using similar inputs as those discussed in our previous work [44]. The power law running term on the right hand side 

is controlled by a linear combination of slope parameters denoted by h. We have checked that with the assigned string 
mass scale, ~ 1 TeV, we can satisfy the above relation by using the indicative values, 6 ~ 20, ^ ~ 5. 

We now turn to the predictions for the chirality conserving contact interactions of four quarks and/or leptons. 
In view of the uncertainties on the renormalization group scale evolution of the gauge coupling constants and the 
coefficient functions, we have chosen to express gg by selecting the electroweak gauge coupling constant, <?(, = 7752, in 
the defining formula Qs = i— M!L^2i!Ii£2 ^ with the proportionality factor set at 77 = 1. 

The information on the branes configurations and on the structure of the associated coefhcient functions is displayed 
in Table I for Cremades et al., [60, 61] model. The volume of cycles may vary substantially from one brane stack to 
the other, so it is important to keep track of the data assigned to the brane configurations DABC which affect the 
normalization factor of the local operators. Since we could not find any analytic approximation that yields reliable 
estimates for the coefficients g^^i in the appropriate range of m^r values, we have numerically evaluated Eq. (11.20) 
by following the same procedure described in our previous work [44]. The contributions from the localization and 
classical partition function factors in the .x-integral arc evaluated by numerical quadrature after removing the massless 
and momentum mode contributions by subtracting the leading terms near the end points x ^ and x ^ 1, according 
to the prescription described schematically by Eq. (11.15). The series simimations over the world sheet instantons 
must be carried out to large enough orders, max(pA) — 6 — 10 and the x-integral must be evaluated with care. (We 
have made use of the Mathematica package.) 

The munerical values of the flavor diagonal {A.F = 0) coefficients, , obtained with the reference set of param- 
eters at the three values of the effective radius parameter, mgr = 1, 2, 3, are displayed in the last three columns of 
Table I. These results all refer to the Ai^ = configurations with coincident intersection points, eA = = (0,0,0). 
The presence of strong cancellations from competing terms appears clearly from the fact that the coefficients change 
sign with variable mgr. The results for the various flavor and chirality configurations are seen to cluster around two 
group of values associated with the pure and mixed chirality configurations, /| ^ and /|/^, which also correspond 
to the cases with equal and unequal angles. The coefficients in the unequal angles group, e|e|j, q\d\, uj^e'^, are 
roughly equal and separated by a gap of about a factor 5—10 from the coefficients in the equal angles group, 
qj^, e%, e|j, q\e\, dj^ej^- Since the exchange of massive vector and axial vector bosons contribute pure and 

mixed chirality amplitudes of same and opposite signs, respectively, wc infer from the comparison of the coefficients 
with same and opposite chiralities that the string excitations are akin to linear combinations of vector and axial vector 
modes. The flavor diagonal coefficients all feature the power law growth with the radius parameter, \gf^fi \ oc (mgr)^/^. 
It is worth noting that the approximate representation in Eq. (II. 9), subsuming the contributions near the boundaries 
of the x-integral, would not reproduce the observed power law dependence of the coefficients on mgr. 

The observed regularities in the coefficients are partly accounted for by the symmetric character of the brane setup 
in the model at hand. Since the branes a and d are always parallel in all three planes, equal interactions are found for 
qj^ and ej^ and for dj^ and e|j. Numerically close values are also found for, u'j^, d'^^ and e^'', due to the fact that the 
brane angles in the various configurations are equal up to permutations of the complex planes. The mixed chirality 
coefficients are larger because they involve brane configurations with two sets of unequal angles. 

For a clearer assessment of the dependence on rUgr, we display in Figures 1, 2 and 3 plots of the coefBcients g^j^i ■ 
We study here the sensitivity on the geometric parameters by varying these one by one with respect to the reference 
set. Note that the longitudinal distances ca, es are associated with flavor change while the transverse distances are 
associated with gauge symmetry breaking. We group the configurations into three classes corresponding to the flavor 
change: AF = : = es = 0: AF = 1 : = 0, es 7^ 0; and AF = 2 : 7^ 0, es ^ 0. Figures 1 and 2 refer to 
the unmixed chirality configurations with equal brane angles and Figure 3 to the mixed chirality configurations with 
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TABLE I: Chirality conserving contact interactions of two quark and/or lepton pairs of fixed chiralities. In Column 1, we specify 
the four fermions configuration; in Column 2, the target space polygon DABC which realizes the world sheet embedding; in 
Column 3, the interbrane angles 9, 9' associated to the two pairs of conjugate fermions f, f; in Column 4, the quartic traces 
over the gauge matrices including the direct and reverse permutations, T1234 ~ T1234 + 74321, multiplied by the Dirac spinor 
matrix element, iSi234, using the conventional ordering for the flavor and color indices i,j,--- and a,/3,--- of the incoming 
fermions, fi,i,a{—ff)f2,j,i3{ff)f3^k,-i{—9')fi,i,5{6'), with the dependence on color indices defined by, Tag-ys = ^{SagS^s +Sag5^g); 
and in Column 5, the numerical predictions in Cremades et ai, [60, 61] model for the flavor diagonal coefficients gfifi = 
{m^gKA/ ^^J^,g\) in Eq. (11.20) at the three values of the compactification scale parameter, mgr = 1, 2, 3. 



Fermions f^f'^ (DABC) 6, 9' ^^1234 '5i234 


9fff'{L,) 


27V ^ 

(eL)^(eL)^ ibdhd) Odb, Odb -jrr{^^Lj^^'2L)(e3Ljae4L) 
(elf {elf {dc'dc') 9,,^, 9,,^ ^(e2«7^eifl)(e4«7Me3fi) 
{eLfielf {dc'db) 9,,d, 9m ^i(e2fl7''eifl)(e4i7^e3i) 


+0.049 - 0.270 - 0.952 
+0.049 - 0.270 - 0.974 
-0.405 - 1.92 - 5.74 


{qLfiqLf ihaba) 9ab, 9ab ^T^TagjsiqiL'y'"q2L){q3Lyij.q4L) 

iuiy{ulf {acac) 9ca, 9ca '^^Tal3js{U2RYuiR)(U4,R'y^U3R) 

{dlfidlf {ac'ac') 9c'a, 9^,^ '^§^To,g^s{d2RYdiR){dAR'yf,d3R) 

{qLf{ulf (abac) 9ba, 9ca -J^^Tacg-fsiq2L'y'^qiL)iU4,R'y,iU3R) 

{iLfidlf {abac') 9ab, ^'c'a -j^^Sag~f5{q2L'y'^qiL){diR'yi^d3R) 


+0.049 - 0.270 - 0.952 
0.0422 - 0.286 - 1.015 
+0.049 - 0.270 - 0.974 
-0.388 - 1.85 - 5.54 
-0.405 - 1.92 - 5.74 


{qLf{eLf {babd) 9ab, 9db ^J^|(5c«;3(giL7^(j2L)(e3i7A,e4L) 
{uVfe-lf {cacd') 9ac, 9d'c ^^^Sag{uiR'y^U2R){e4,R'y^e3R) 
{dlf{elf [cacd) 9ac', 9dc' ^^|Ta/3(difl7^d2fl)(e3fl7Me4fl) 


0.049 - 0.270 - 0.974 
-0.165 - 1.43 - 4.15 
-0.0157 - 0.835 - 2.22 



unequal branc angles. The presence of cusp discontinuities in certain curves is due to our use of logarithmic plots for 
the absolute values of the coefficients aimed at representing quantitatively the size of the suppression. 

We see on panel (a) of Fig. 1 that the predictions are spread by an approximate factor 2 — 3 for reasonably 
restricted variations of the shape parameters. Using tilted 2-d tori, or increasing the complex structure parameter 
X^, results in enhanced coefficients, while decreasing results in reduced coefficients. The A_F — 1 coefficients with 
finite eB in panel (6) are suppressed by order 10^^ while the A_F = 2 coefficients with finite tA and es in panel (c) 
are suppressed by factors of order 10~^ — lO"**. The specific dependence on rUsV is a result of the tension between 
the power growth from the overall normalization factor gg oc (mgr)^ and the exponential suppression from Zd- That 
the suppression effect is controlled by the classical factor is clear from the fact that the coefficients have comparable 
values near mgr = 1. 

The plots in Fig. 2 again confirm that < 1 and > 1 lead to reduced and enhanced coefficients. The nearly 
one order of magnitude suppression of the AF = 1 coefficients is independent of ■ That the suppression is weaker 
than expected is explained by the specific feature in the present model that only a single component of the vectors 
eg are finite. The comparable predictions found for cb = (0, 1/3,0) and (0,2/3,0) are explained by the torus lattice 
periodicity. The cancellation effects from the oscillating factors e^*'^^'^ explain both the change of sign from positive 
to negative coefficients and the smooth variation with mgr. 
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The plots in Fig. 3 show that the coefRcicnts in the unequal angle cases are systematically larger than those found 
for equal angles. The dependence of the wrapped cycle volumes on the torus shape parameters spreads the coefficients 
by a factor of 2 — 3. 

Two general features of the predictions are the rapid power law increase with rUgr of the flavor conserving AF = 
coefficients and the hierarchies of order 10~^ and 10"'* separating these from the flavor changing AF = 1 and AF = 2 
coefficients which vary more slowly over the allowed interval for m^r. While the variation of the coefficients with 
msT is not apparent on the result in Eq. (II. 9), obtained by restricting the x-integral to the end point contributions, 
it appears possible to use this approximate formula in order to explain the dependence on the distance parameters, 
eA, es- Examination of the combined contributions from the string momentum and winding modes to the coeflicients 



HH' _ \^A\ V- V- 11/ "/ ^ - ^ ^TT o2^ 



shows that for small finite the leading contribution to the ratio of AF = 1 to AF = amplitudes is of form 
^-PAiAlCAlimsr) ^ t^YiHq tiie = 2 amplitudes include the additional suppression from the oscillating factors, e^'^**^. 

It is interesting to compare our predictions for the contact interactions of four fermions with those obtained in the 
Dp/D[p + 4)-brane models [33]. (The formalism is briefly reviewed in Appendix B.) For the coupling of four modes, 
1(3, 7j)p|(3, 7j)|^, the comparison at fixed 4-d gauge coupling constant of our estimate, «Si234'7i234 ~ 27r(0.05 — 

0.5)(msr)^5i2347i234, with the result found by Antoniadis et al., [33] in the large radius limit, 7i234[0.12 ^1234 + 
0.33 1S1234], reveals an order of magnitude concordance. 

Finally, we compare our predictions with the contributions from the momentum modes evaluated by means of 
Eq. (11.16) for n = 3. The resulting rough estimate, {gfjf )kk — i-K(n-2) ■^»('^ji^)"l^-B| — {'m.s'rY' Cb, indicates that 
the contributions from the string momentum modes are significantly larger than those from winding modes. We 
should remember, however, that the present estimate must be regarded as an upper bound since it relies on the large 
r limit and ignores the form factor suppression. 



III. INDIRECT HIGH ENERGY COLLIDER TESTS 



We discuss in the present section the collider physics applications based on the formalism presented in Subsection II A 
and in Appendix B for the orientifold model of Cremades et al., [60, 61]. Since the distinction between the mass and 
gauge bases is not essential for these observables, all the results in this section are obtained by setting the flavor 
mixing matrices to unity, = 1. 



A. Bounds on contact interactions mass scales 



It is important to distinguish the mass scale A associated to the 2? = 6 operators from the mass scale Mh associated 

to the I? = 8 operators in the 4-d effective Lagrange density quadratic in the energy- momentum tensor [78] , Leff = 
i-^Tni^T^^" . The analyses of available high energy collider experimental data using field theories in extra space 

dimensions are sensitive to values of these mass scales, Mh = 1.5 TeV [78], A = 2 — 6 TeV [79] and A = 1 — 8 TeV [80]. 
In the single _Dp-brane models, the quantum gravity mass scale Mh was found to be parametrically larger than the 
string scale [25], ~ vJtt ' thus making the detection of new physics effects harder. To pursue the comparison 
with the intersecting brane models, it is convenient to consider in place of Mh the closely related gravitational mass 
scale Mc defined in the case with n flat extra dimensions by [11], Mg+V" = ^ = {4:TtGn)~'^- The mass scale Mq 
is related to the fundamental string parameters of single Z)3A-brane models as [25], (MG/nis)^ = 167T/g\. Repeating 
the same calculations for intersecting F)6A-branes gives us the modified formula in the large radius limit 

.Mg 8 ^ ^Stt 

^ms' g\{m.rf\m'^-n'^U'\/Ui 

The strong dependence on the geometric parameters indicates the interesting possibility that Mq may assume lower 
values in multiple brane models. 

We now discuss the constraints on the string mass parameter, rris, inferred by comparing our predictions for the 
chirality conserving contact interactions of 2? = 6 with a subset of the available experimental limits [81]. For the lepton 
and lepton-quark configurations, e\, e\q\ and the quark configuration, g^, respectively, we evaluate the bounds on 



(III.1) 
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the string scale parameter, = gfiii^^^li^ + 5^f)g\/ Ka^ at fixed TOsT, by setting the gauge coupling constants, 
g\/KA at 5I/2 ~ 0.213 and g1 ~ 1.44, respectively. Using the numerical values of the coefficients in Table I, we 
obtain the following bounds on mg for the choice of three representative experimental limits on the mass scales: 

Aj^(eeee) > [4.7, 6.1] =^ rus > [0.69, 2.1, 3.9] 
Aii(eegg) > [23.3, 12.5] => > [2.4, 3.0, 5.6] 

Afita) > [2.7] ^ ms> [0.72, 1.7, 3.2], (III.2) 

where all masses are expressed in TeV units and the three entries refer to the values m^r = [1, 2, 3]. 

We next consider the constraint from the enhanced supernova cooling through the reaction producing right handed 
neutrino-antineutrino pairs by quark-antiquark pairs, q + q ^ + P", which is allowed as long as the contributions 
to the neutrino Dirac mass arc bounded by the supernova temperature, m,^ < Tsn — 50 McV. The lower bound on 
the mass scale in the effective Lagrangian, Leff = ^gl^a (g7M75g)(^fi7^^fl)> is found for the SN1987A to cover the 
range [82], A'^"^ > (90 - 250) TeV. For concreteness, we set our choice on the lower bound, A'''^ > 200 TeV. Using 
the numerical predictions in Table I and assuming the approximate equalities between the chirality basis amplitudes, 
qj^i^]^ — qj^uji and dj^iyj^ ~ dj^, we obtain the numerical estimate for the coefficient, [gq^y^^ ~^gq,^ii) — [0-66, 1.26, 2.14] 
for rrisr = [1, 2, 3]. The resulting bounds on the string mass scale road, nis > [43, 82, 139] TeV. For reference, 
we note that the comparison with the contribution from the string momentum modes yields [50] the weaker bound, 
> (5 - 10) TeV. 



B. Bhabha scattering cross section 

Useful constraints on the new physics are set by the experimental data at the high-energy colliders involving the 
two-body processes [84] of Bhabha, Moller and photon pair scattering and fermion-antifermion pair production. The 
absence of significant deviations from the Standard Model predictions has led the statistical analyses of data to set 
exclusion limits on the free parameters. The global fits to the combined high-energy collider data based on the single 
£>-brane model, with the gauge factors treated as free parameters, yields [83] rUg > (0.69 — 1.96) TeV. 

Wc focus here on the Bhabha scattering differential cross section for which high precision mcasTircments along with 
higher order calculations of the pertubation theory corrections are due in the future. Experimental data has been 
collected by the LEP collaborations [85-87] . The studies based on single Dp-brane models, describing the ratio of the 
string to Standard Model differential cross sections by the approximate formula [25] 

5^ 

yield by comparison with the experimental data at the center of mass energy ^/s = 183 GeV the 95% confidence 
level exclusion limit on the string mass scale, rris > 410 GeV. Similar bounds are found in the analysis [88] including 
the experimental data at ^/s ~ 188.7 GeV. At the higher energy, y/s = 1 TeV, the 95% confidence exclusion limit 
obtained under similar conditions should extend the sensitivity reach to [25], nis > 3.1 TeV. 

We now present our predictions for the Bhabha scattering differential cross section evaluated with Cremades et 
al., [60, 61] model by adding the contributions from the contact interactions in Eq. (11.15) to the Standard Model 
amplitudes, using the formalism detailed in Appendix B 2. In Fig. 4 and Fig. 5, we show plots of the ratio R{cos9) as 
a function of the scattering angle variable, cos^, for the center of mass energies, y/s = 183 GeV and y/s = 500 GeV, 
respectively. The selected set of values for rus are different for these two cases, as dictated by the fact that the string 
corrections scale as s/mj.. 

For a qualitative comparison with experimental data, we note that the LEP data points for the ratio at -^i = 0.183 
TeV are spread over the interval of cos^ e [—1, -|-1] inside the band limited by the horizontal lines at, i?(cos^) = 
(1.0 ± 0.4). As for the single D-brane model prediction [25] in Eq.(III.3), this is represented by a nearly straight line 
which slopes from 1.06 to 1.0 as cos 9 increases from —1 to 1. The exchange of string Regge and winding modes are 
seen to give a small reduction of _R(cos(?) near the forward scattering angles, cos^^ ^ 1, gradually turning into a large 
enhancement near the backward angles, cos^ — 1. This implies a change from a relative negative sign to a positive 
sign at some intermediate angle in the interval cos^ G [—1, -1-1]. The contributions grow rapidly with increasing 
msT. Requiring the predicted ratios in Figs. 4 and 5 to remain bounded inside the interval R{cos6) G [0.8.1.2] 
for cosO e [—1, -|-1] imposes lower bounds on the string scale which cover the ranges, rUs > (0.5 — 3.) TeV and 
ms > (2. — 5.) TeV, respectively, for the interval of values mgr G [1,3]. 
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We have also performed a more realistic calculation of the Bhabha scattering differential cross section in which the 
total regularized string amplitudes are obtained by subtracting by hand the pole terms from exchange of massless and 
massive momentum modes for the j + Z gauge bosons, while adding up the contributions from the physical j + Z 
pole terms, using the prescription in Eq. (11.15). To ease the numerical calculations we have only evaluated the pure 
chirality amplitudes LL, RR, while assuming the mixed chirality amplitudes LR, RL to be proportional to these. 
The model dependence on the ratio of pure to mixed chirality amplitudes thus resides in the adjustable proportionality 
constant, x = Ge^et/Gej^et — G^^ /G^^ = G^^/G^^, which we have taken to vary inside the interval, x = 5]. 

The ratio i?(cos 9) of the predicted differential cross section to that of the Standard Model is plotted in Fig. 6 at 
the center of mass energies, -y/i = 183 and 500 GeV (left and right hand panels) for the two values of the string scale, 
nis = I, 2 TeV and mg = 2, 4 TeV, respectively. The comparison of the curves /, //, and similarly of the curves 
III, IV, measures the sensitivity of R{cos6) with respect to the string scale rUg. On the other hand, the comparison 
of the curves /, ///, and similarly of the curves //, IV, measures the sensitivity with respect to the mixed chirality 
amplitudes. The large spread of predictions with variable x and nis shows that Bhabha scattering can usefully test 
the model dependence. The results from the present complete calculation agree qualitatively with those in Figs. 4 
and 5, although the change of slope and subsequent rise of the ratio with decreasing cos 6' are generally less steep. We 
conclude that the representation of string amplitudes by contact interactions is reliable for the considered incident 
energies. 

IV. FLAVOR CHANGING NEUTRAL CURRENT PROCESSES 
A. Direct and indirect flavor changing eflfects 

The flavor changing neutral current observables are determined by the non-diagonal elements of the mass basis 
coefficients of contact interactions, Gij.ki- These receive direct flavor changing contributions from the four point 
string amplitudes and indirect contributions from the linear transformations linking the gauge and mass bases of 
the fermions. Without further input information on the flavor structure, it appears impossible to infer quantitative 
constraints from a comparison with the flavor changing observables. 

A natural description of the fermions flavor quantum numbers is provided by the basis labeled by the branes 
intersection points. The geometric constraints on string amplitudes, as stated by Eq. (11.23), directly translate as 
selection rules on the flavor amplitudes in this basis. For Cremades et al., [60, 61] model, however, these rules turn out 
to be trivial ones, owing to the fact that the intersection points for the fixed chirality modes lie at finite distances apart 
only in single complex planes. Since the conditions involve at least one shift vector defined modulo 1, no zero entries 
are enforced either on the trilinear Yukawa interactions, XjjfifjH, or on the four fermions interactions, Gij^klfifjfkfl- 
This property is also responsible for the separable structure of the Yukawa coupling constants, Xjj = aibj, implying 
that the mass matrices of quarks and leptons all have unit rank. Although the fiavor non-diagonal coefficients are 
generally finite, the AF = 1,2 operators associated with the configurations, i = j ^ k = I and i ^ j ^ k, are strongly 
suppressed with increasing mgr by the classical partition function factor for longitudinal distance parameters, sa, (b 
of 0(1). However, the fact that the tree level string amplitudes depend on the relative distances between intersection 
points, efj and e^, introduces certain restrictions on the flavor structure of the coefficients Gij^kl- 

We here focus on the two flavor changing neutral current observables associated to the mass splitting of CP 
conjugate pairs of neutral mesons made of quark-antiquark pairs, P = QiQj, P = QiQj, and the three-body decays 
of leptons, Ij li + Ik + h- For simplicity, we assume that the Standard Model contributions to these observables 
are negligible relative to those from the contact interactions, so that we can directly use the experimental limit to 
derive bounds on the string scale. Convenient formulas for the contributions to these observables from the chirality 
conserving local operators have been obtained in [89] for models with extra U{1) gauge symmetries. We follow closely 
the formalism developed in the latter work, while accounting for the fact that the dependence on the color quantum 
numbers is different in our case. The observables for the real and imaginary parts of the P° — P° mass splitting. 
Amp = -3fJ(< Po|Lg|=2|pO >) and ep = -9(< Po|Lg|^^|P° >)/(2^/2mp), are given in our notations by the 
explicit formulas 

Amp = -2mpFl,[l^{Gf/:,. + G^%) - a^K(G^.f^,.)]. 

where a'- - ~ h + |( — ^ — Note that the definition for the indirect CP violation observable applies specifically 
to the K — K system, ep = ck- The quarks fiavor indices i, j are set in accordance with the conventional assignments 
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for the neutral mesons, K^{sd), B^ibd), B^ihs), D^{uc), with Fp denoting the mesons two-body leptonic decay 
couphng constants. We have evaluated the hadronic matrix elements of the four fermion local operators for the 
pseudoscalar mesons to vacuum transition by making use of the vacuum insertion approximation. The bilinear axial 
current hadronic matrix elements are determined through the PC AC hypothesis in terms of the measured parameters 
Fp. Using the conventional definitions for the matrix elements of quark bilinear current operators 

< Q\SocLl^d0L\K''{p) >=-< K''{p)\-S^Ll^df)L\(i >= i ^J^ P^Sap, 

< 0|s„75rf/3|i^°(p) >=< ^0(p)|s,75rf/3|0 >= -i-—!^!:^ -6,0, (IV.2) 



one can write the matrix elements of the relevant quadratic operators as 



Qal3-/5-^al3-fS - ^ , Qaf3j5-''Oc0'r5 ^ [1 + ^[ ) \, (IV.^j 



where the saturation of color indices displayed in the last entry above uses the tensor, Tap-ys — \{5afi5yS + SaS^-yis)- 
The factor a-^ accompanying the coefficient Gij%, differs from that quoted in [89], aij = ^ + ™^ which 

refers to the dependence on color indices involving the diagonal tensor, Sa^S^s- Similar formulas hold for the B and 

D mesons. 

The contributions from contact interactions to the lepton number violating three-body decay rates of charged 
leptons are given by [89] 

Tn 

T{ej - e, + e, + ei) = ^[i\Gf,%\' + \Gf,%\' + \Gf,%,\' + \Gf,%,\') + (L - R)]. (IV.4) 

The partial rates for the pair of decay reactions, fi~ — > e~ + e+ + e~ and r~ — > e~ + e+ + e~ , with j = 2, 3 and 
i = k = I = 1, are described by the simplified formula 

r(e, ^ e. + e. + e.) = ^^i^f[2\gf/.f + 2\g^f + \g^f + Ig^f], (IV.5) 

where we have included an extra symmetry factor i in order to account for the pair of identical charged leptons in 
the final states. 

The bounds on the string scale implied by the mesons mass shifts and the charged leptons decay rates are expressed 
by the explicit formulas 



rus > 



^ ,2'!Tg\^i 



777^ 

S84Tr^r{ej ^ a + a + Ciy \ 



1/4 



(IV.6) 



In view of the partial information that we dispose on the matrices Vfj and the complicated summations over the 
flavor basis amplitudes, we choose to perform an approximate calculation motivated by the specific flavor structure 
of contact interactions for the model at hand. We only retain the coefficients Qij^i with i = j, k = I, denoted by 

gfl^ = g^^i^, and neglect the distinction between intersection points, by assuming the diagonal terms gu to be 
independent of i and the non-diagonal terms to be symmetric, gij = gji. Applying now the unitarity conditions on the 
quarks transformation matrices, V^, allows us to express the mass basis coefficients for the neutral meson observables 
in the form 

~9^jfi - (^f^^Ai + ^'^iKj + hrAj){gf,"' - gg"') + 2hifh,^{gg"' - g%"') + (K^^ + hijh.j){gfj"' - gf^"'), 
[hij = V^.jV^l, = V^lV^.., i^j^k] (IV.7) 
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where h'^j, h'^j arc given by same formulas as hij, hij with H ^ H' . Assuming further that the non-diagonal elements 
are independent of the specific pair, z, j, so that 1712 = 513 = (7237 leads to the factorized form 

[V^ = 2hijh,j + 2h,fh,j + hijh'i^ + h^jh,j\. (IV.8) 

We have considered the alternative approximation defined by assuming that the diagonal coefficients ga = are 
independent of i and the non-diagonal ones satisfy git = gki — gnd, but without imposing the unitarity conditions. 
The resulting form of the AF = 2 coefficients reads 

k kjtl kjtl 

For the coefficients gf^ff entering the charged Icptons decay widths, we use the same assumptions as in the calcu- 
lation done just above to obtain the simplified formula for the mass basis coefficients 

= E ^HjkVh*ikVwikVh*ik, Vji = E ^HjkVh^ikVljHiVh^i]- (iv.io) 

k kjtl 

B. Results and discussion 

We have numerically calculated the mesons mass shifts and the three-body leptonic decay rates for Cremades et 
al., [60, 61] model with the reference set of parameters described previously and the longitudinal distances set at, 
£a = (0, 0, 0), es = (0, 1/3, 0). Somewhat arbitrarily, we choose to set the various flavor mixing matrices equal to the 
CKM matrix, V^^ = Vckm, [H = L,R]. The following input data, expressed in GeV units, are needed to calculate 
the mass shifts: 

m„ = 3 X 10"^, TOd = 7 X 10"^, rus = 0.095, TOc = 1.25, rub = 4.20; 
Fk = 0.1598, rriK = 0.497, Atuk = 3.483 x 10'^^, = 0.0979; 

Fb = 0.176, ruB = 5.2794, Ams = 3.337 x IQ'^^, V^" = 0.000125; 

Fb, = 0.176, niB, = 5.367, Atub, = 1.145 x 10"", V^'' = 0.00347; 

Fd = 0.2226, ruD = 1.8645, Amo = 4.607 x 10"", = 0.0979. (IV.H) 

The following input data, expressed in GeV units, are needed to calculate the charged leptons three-body decays: 

rUe = 0.511 X 10-^ = 0.1056, r(;U ^ e -|- e -|- e) < 3.29 x W'^^ 5RV21 = -0.198633, KV21 = 0.198722; 
rrir = 1.777, r(r ^ e + e + e) < 4.529 x 10-^^ ^{Vsi) = 0.00641, ^{V's-^) = -0.00660. (IV.12) 

The results obtained with the flavor mixing described by the approximate formulas in Eqs. (IV.8) and (IV.IO) are 
plotted in Fig. 7. The use of Eq. (IV. 9) gives similar results. We see that the bounds on m,,, at fixed m^r, increase 
with increasing mgr according to the approximate power law, (mgr)^/"*. Wide disparities appear between different 
cases mainly because of the flavor mixing factor. The most constraining observable, corresponding to the K — K 
mass shift, yields the bound, rris > O(IO^) TeV. Relaxing our assumption that the Standard Model contributions are 
negligible can only strengthen the bounds on m^. However, one may expect significantly weaker bounds if the flavor 
and mass bases were not too strongly misaligned so that the flavor change is dominated by the direct contributions. 
For instance, using the order of magnitude predictions for the off-diagonal coefficient gfj^i in panel (c) of Fig. 1, with 

= 1, would reduce the bound on rus from the K — K mass shift by a factor of order 10~^ — 10~^. A careful 
treatment of the flavor structure of the model would be needed in order to make a more definite statement. 

For comparison, we note that the bound from the K — K mass shift obtained in [49, 50] by using the approximate 
representation of Eq. (II. 9) for the string momentum modes reads, mg > 100 TeV. A similar gap exists for the other 
fiavor changing observables. However, these results were obtained by setting, m^r ~ 20, which lies well above the 
allowed interval. 

We comment briefly on the CP violating observable, ck, which is set by the experimental data for the n + n 

decays to the value, jex je^p ~ jryooj = (2.285 ± 0.019) 10~^. Since the coefficients g^^i are real, the prediction for 
eK depends in a crucial way on the flavor mixing matrices. In our treatment of the indirect flavor mixing leading to 
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Eq. (IV. 8), the predictions for and the mass shift scale as, |e/f/Am/f | = |3(V^)/5R(V^)|. Since the CP violation 
effects enter the CKM matrix through the second and third fermion generations, is real and hence = 0. The 
alternative prescription for the flavor mixing described by Eq. (IV.9), with the matrices still identified with the 
CKM matrix, yields an uninteresting small bound on m^. 

It is also instructive to compare with the split fermion models. The bound from the K — K mass shift [56], 
Mc > /3[V^F(/9a)]i ~ (100 - 600) TeV with /3 = 1125 TeV and F{pa) ~ (2. -8.), is of same magnitude as ours, while 
sampling over the parameters which control the indirect flavor mixing effects give the ability to suppress the bound by 
factors of 10 — 100. The description of exchange contributions in these models differs from that in intersecting brane 
models where the flavor hierarchies originate in the instanton contributions, rather than the wave function overlaps, 
and the parameter space is more restricted. In addition to the extra dimension size parameter, Mc = 1/r, the split 
fermion models [56] introduce the scaled localization width and fermion separation parameters, p = a/r and a = ^yja, 
which qualitatively identify with the string theory parameters, p ~ Mc/nis and a ^ erus/Mc = emgr, assuming 
a ^ l/rris. We conclude from this indirect comparison that the wide hierarchies, a £ [0. 15] and p £ [10~^, 10"**], 
which are needed in split fermion models to weaken the bounds on are not favored by the analogous intersecting 
brane models. 

Finally, we present the result of an indicative study of the tau-lepton hadronic and semi-hadronic decays, t n + p 
and r — > TT + TT + /U, based on the analysis [90] of the effective interaction for the associated subprocesses, ALeff = 

■]^y{pi^T){qjf^q) + H.C., which yields the bound, A'^'^ > 12 TeV. Using the predictions in panel (c) of Fig. 1, we 

deduce the bound nis = ^'^'^Qaigu^)^ > 0.5 TeV at rngf = 1. At larger m^r, neglecting the flavor mixing effects leads 
to useless small bounds due to the strong suppression of the flavor non-diagonal string amplitude. 



V. SUMMARY AND CONCLUSIONS 



We have discussed in this work collider and flavor physics tests of the four fermion tree level string amplitudes in 
intersecting brane models. Although the study was specialized to the isolated orientifold premodel of Cremades et 
al., [60, 61] realizing the Standard Model, this is a good representative of the families of string models selected in 
current explorations of the landscape of open string vacua. Based on a qualitative examination of the predictions for 
the gauge coupling constants, we also verified that it is compatible with a TeV string mass scale. 

The string theory predictions depend on two free mass parameters, rUs and M,-; = l/r, along with the geometric 
shape parameters of the internal torus and known inputs for the electroweak gauge bosons masses and gauge 
coupling constants. The necessary condition for weakly coupled open strings imposes the restricted variation interval, 
nigr = rris/Mc G [1,4]. We have studied the four fermion contact interactions from exchange of string Regge and 
winding modes in various configurations of the quarks and leptons, paying special attention to the gauge group 
structure and the contributions from world sheet instantons. The general features of predictions for the contact 
interactions may be briefly summarized as follows. The size of coefficients present regularities which reflect in part 
the symmetric configm-ation of the brane setup. The sensitivity of predictions to the tori shape parameters leads to a 
moderate sensitivity of the flavor diagonal coefficients AF = on geometric parameters which spreads predictions by 
a factor 2 — 3. The widest disparities occur between the mixed and unmixed chirality amplitudes. Two characteristic 
features reside in the strong growth of the flavor diagonal coefficients, (wsr)^/^, and the strong suppression of the 
flavor non-diagonal AF = 1, 2 relative to AF = by factors of order 10~^ and 10^'', due to the classical partition 
function factor when the distances between intersection points relative to the wrapped cycles radius are of e = 0(l/3). 

The Bhabha scattering differential cross section is an important high precision observable for which the theoretical 
and experimental uncertainties are expected to reach O(10~'^) in the future. We have considered a qualitative 
comparison with the LEP data which leads to bounds on the string mass scale of TeV order. These are expectedly 
stronger than the bounds obtained in the single brane model [25] where the local operators have dimension 8. It 
should be useful to pursue a systematic study for the set of 2 — 2 body processes including the DrcU-Yan lepton pair 
production and the parton subprocesses with initial states, for e + q, q-\- ^ and q + (( . 

We have also considered the direct contributions to the four fermion contact interactions from string Regge and 
winding modes to a subset of flavor changing neutral current observables using an approximate description of the 
indirect flavor mixing effects where the direct and indirect flavor changing effects factorize. The K — K mass splitting 
yield the strongest constraint, TOg > 10'^ TeV. This bound, as well as other ones deduced from flavor changing 
observables, are an order of magnitude stronger than those obtained from the contributions to contact interactions 
due to the string momentum modes [49, 50]. It is fair to say, however, that this conclusion is at best qualitative 
since the two calculations rely on different inputs and approximations. To obtain more realistic estimates, the highest 
priority should be set on obtaining realistic inputs for the flavor mixing matrices which match the predictions for 
the fermions mass matrices to observations while improving on the restrictive rank 1 property of the fermions mass 
matrices in the model at hand. 
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APPENDIX A: BRIEF REVIEW OF OPEN STRING SECTORS IN TOROIDAL ORIENTIFOLDS 

We consider type // string theory compactified on factorisable toroidal orientifolds, T^/CITZ, [T^ = 0/=! where 
the involution symmetry, TZ = 11/=! '^i^ ^'^^^ '^^ orthogonal basis of complex coordinates, = {X( + jX|)/\/2 
in the Tf, [I = 1, 2, 3] complex planes as reflections across the real axes, TZi ■ = X^ . We restrict to the subset of 
factorisable three-cycles in the integer homology vector space, 11^ e H^{T^, Z), represented in terms of the homology 
of one-cycles with lattice dual bases, [a^], [b^] E H^(Tf, Z), by the three pairs of integer quantized winding numbers, 
(n^, m^). These three-cycles are represented in the orthogonal coordinate system of the three complex planes of 

by 

[n^] = K)' K = K- K^l^ V' ^ + iu^ = -4 = -bi + ^,xi = 4] 

where denote the Tf tori complex structure moduli whose real parts are subject to the restriction, 6/ = 0, i, 
for orthogonal and tilted tori, respectively. We have denoted by r{, rl the radius parameters of the one-cycles of 
Tj projected on the pair of orthogonal axes. For non-orthogonal 2-d tori, , we choose to work with the case of 
upwards tilted tori, where r| refers to the one-cycle along the imaginary (vertical) axis of the complex plane and r{ 
to the projection of the dual one-cycle radius along the real (horizontal) axis. The orientifold 06-planes arc the loci 
of points fixed under TZ which extend along the three uncompactified dimensions of Minkowski space-time, M4, and 
wrap the three one-cycles, (n^,m^) — (1,0). Both the 06-planes and Z)6-branes are sources for the closed string 
RR modes seven-form, C7, with RR charges determined by the winding numbers of the wrapped three-cycles. The 
divergent tadpoles of RR modes due to the 06-planes in the one-loop closed string (Klein bottle surface) amplitude are 
assumed to cancel against the tadpoles in the open string (cylinder and Mobius strip surface) amplitudes contributed 
by introducing K parallel stacks of branes -D6^, [/i = 1, 2, • • • , if = a, 6, • • •]. To the D6^-brane stack wrapped 
around [11^] = (n^, m^, is associated the orientifold mirror image DSju'-brane stack, wrapped around the image cycle 
[n^,] = {n^n' , rnj^r) = (n^, — ?7i^). For toroidal orientifolds, the RR tadpole cancellation conditions are of form 

K K 

N^^ny^n^ + -Oop = 0, ^ N^^m'^m^ = 0, [Qop = T'^P'^U = T32] (A.2) 

where the summations run over the orientifold equivalence classes, counting mirror pairs as single units; I, J, K € 
[1, 2, 3] run over the distinct permutations of the complex planes indices; fp = 2^~p denotes the Op-planes multiplicity; 
and the upper and lower signs of the orientifold charge Qop refer to the SO and Sp (orthogonal and symplectic group) 
orientifold projections. 

The Z)6a-brane location in the Tf complex planes is described by an oriented vector tilted relative to the 06- 
plane (along the real axes) by the angles, ttO^ = arctan(TO^/n^[/2 ). The Z)6-branes serve as boundaries for the 
end points of open strings which carry their perturbativc excitations. The open string sectors, (a, 6) and (a, 6'), 
associated to the two pairs of branes, D&a/D&b and D6a/D6b', are assigned the interbrane angles, 9ab = Oh — Oa and 
9 ah' = 9}yi —Oa^ —6}y — 0a- Wc usc notational conventions where the brane-orientifold and interbrane angles vary inside 
the intervals, O^i, & [—1,-1-1] and 9^^ G [0, -|-1], [/ = 1,2,3] with positive sign angles associated to counterclockwise 
rotations. Transforming back to values of the angles inside these ranges requires geometric information on the signs of 
winding numbers. The brane pairs a, b intersect at fixed numbers of points determined by the topological invariants 

lab = Yliniml - mini), I^'t = Uini^ml - m[.nl) = J\(nifhl + mini). (A.3) 

The low-energy dynamics on a single isolated stack of A'^ Z)6-branes in the 4-d space-time is approximately 
that of a gauge field theory with gauge group U{N), supersymmetry A/" = 4, and a certain content of massless 
modes associated with the branes moduli. The open string sectors for the D6a/D6b-hrane pair supporting the gauge 
symmetry U{Na) x U{Ni,) include: (1) The diagonal modes, {n,iJ.), [/i = a, b] which carry the adjoint representations; 
(2) The orientifold twisted modes, {fj,,ijf), which carry the symmetric and antisymmetric representations A, S of 
U{Nfj,) with the multiplicities, ^(/^^' ± Ifioe), [Iiioe = [tt^] • [t^oq] = n7(~™i")]; ^■^d (3) The non-diagonal 'twisted' 
modes, (a, b) ~ (6, a)^ and (a', b) ^ (&, a')\ which carry the bifundamental representations, Iab{Na, Nb)(B Ia'b{Na, Nb). 
The equivalence relations between open string sector sectors, (a, b) ~ (6', a') ~ (6, a)^, (a, b') ^ (6, a') ^ {b' , a)^ , where 
the dagger stands for the complex conjugation of the states space-time and internal group quantum numbers, lead to 
interpret the intersection numbers lab or la'b of negative signs as multiplicities for the modes with conjugate chirality 
and group representation, \Iab\{Na,Nb) or \Ia'b\{Na, Nb). 
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The key condition to preserve Af = 1 supcrsymmctry on the branes 4-d intersection is that the wrapped three- 
cycles be of special Lagrangian type. For Calabi-Yau manifolds, these are the cycles whose real number valued volume 
integrals are calibrated by the holomorphic three-form, e^'Tis, characterized by a fixed choice of the angle parameter, 
ip. In close analogy with the conditions selecting in closed string theories the holonomy of subgroups of the internal 
space manifold symmetry group SO{Q), one preserves jV = 1 or jV = 2 supersymmetry to the extent that the rotation 
matrices relating the branes to the orientifold planes belong to the group SU{3) or SU{2) [4]. The requirement that 
the wrapped cycles Ha and 11;, are calibrated by the same three-form so that the £)6a/-D6b-brane pair preserves Af = 1 
supersymmetry, amounts to the conditions, 0^ ^ = mod 2. In terms of the spinor weights r^^^ of SO{8) for the 
16 supercharges conserved in the bulk, J\f = 1, 2 supersymmetry arises when a single or a pair of spinor weights r(^a) 
solves the equations, Y2t=i ^"aj^^J = 0; = 1; ' ' ' i 4; ^ = 1, ■ ■ ■ , K] where the intersection angle in M4 is set here to 
zero, 6^ = 0. In the basis of independent spinor charges for SO{8) defined by 

r(i) = (- + +-), r(2) = (+ - +-), r(3) = (+ + — ), r(4) = ( ), (A.4) 

with the notational convention, r°^^ = {— + -\ — ) = (— ^, ^, ^, — [a = 1, 2, 3, 4] the three special angle configurations, 
(0, 6a, ±^a), (^b, 0, ±0b), {Oc, ±^C) 0), preserve the = 2 supersymmetries associated to the pairs of spinor weights, 

(•'•(2),r(3)\ (r(i),r(3^\ /r(i),r(2)\ 
Vr(i),r(4)/" Vr(2),r(4)/' \r (3-, ,r J ' 

We discuss next the dependence of intersection points on geometrical data for the angles and transverse separations 
of DG^/DGs-brane pairs that do not necessarily pass through a common point of T^. Suppressing the index I of 
the 2-d tori Tf, for convenience, we parameterize the DG^-brane wrapped around the one-cycle of with winding 
numbers {nA,mA) by the equation 

Xa{^a) = {La^a + Qa+Pat + dAtA)ei, (A.5) 

where 

La = «A + ruAT = nA+ imAT2, Ia = -toaT2 -|- ifiA, nA=nA+ wati, 

[t = Ti+iT2 = — , ei = 27rri, 62 = 27rr2e*"] (A.6) 
ei 

with T denoting the torus complex structure modulus, dA the transverse displacement from the origin, qa, pa & Z 
the lattice displacements along the basis of dual cycles, and ^a, dA ^ R parameterize points along longitudinal and 
transverse directions. We have formulated the problem here in the case of sideways tilted torus, which is related to 
the case of upwards tilted torus considered in Eq. (A.l) by the transformation, r = —1/U. (The formulas for the 
upwards and sideways tilted tori are also related by the substitutions, n ^ m, rn — » — n, ei — > 62, 62 ^ — ei.) We 
now introduce a similar equation for D6b ■ XB{S,B)/ei = LbS,b + Qb + Pbt + dstB, and require the condition, 
-^a(^a) = -^b(^b)- The resulting pair of linear equations for the real variables ^a and ^s, 

( "A -fiB \ f Ca \ _ _ f Qab \ / A 

[niA -TUB J Ub)^ [PabJ' ^^-^^ 

where 



Qab =Qa-Qb, Pab = Pa- Pb, Qij. = qij.- d^m^T2, Pij,=Pij. + d^h^/T2, [11 = A, B] (A.8) 
is solved in matrix notation by 

f Ci \ ( TUB -riB \ ( Qab \ . ^ - ^abLa y _ k'sA^B 

K^bJ IabK^a -nAjyPABj ^ Iab ' Iab ' 

[Iab = nAmB - niAUB, kAB = {fiBPAB - it^bQab) mod(/AB), k'g^ = {uaPab - mAQAB) mod(/AB)](A.9) 

The Iab solutions for ^a, are in one-to-one correspondence with the pairs of integers kAB, ^ba- The shift vectors, 
Wba = Xb — Xa = {k'gj^LB — kABLA)/lAB, which link the positions of a given intersection point along the pair of 
intersecting branes {B,A) in the complex plane of T^, belong to the lattice coset, A^a/A, where Aba denotes the 
grand torus lattice generated by the cycles Lb, La, and A the torus lattice [67]. 

We next discuss the vector space of open string states, {a,b), [a = b, a ^ b]. The state vectors, 
I (a, b)k, N, r(„) , (A, ij) > a|^^^ , are described by the four momentum fc; the coordinates oscillator number N; the weight 
vector r(Q-) of the S0{8) group Cartan torus lattice; and the gauge group multiplct component A. The CP gauge fac- 
tors, Xa^j ' = 1) ■ • • ) A^a! i = 1) • • • ) ^b] are matrices whose array and column indices i = (1, • • • , Na), j = (1, • • • , Nb) 
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label the coincident branes inside the stacks a, b oi size Na, Nf,. In orientifolds, the combined systems of mirror brane 
pairs, (a + a', 6 + are described by single matrices, A^_^^-, [i, j G (a + a', 6 + 6')] of size (2Na + 2Nb) x (2Na + 2Nb). 
These are conveniently represented by 2 x 2 block matrices with sub-blocks labeled by a, a', 6, 6'. The diagonal open 
string states (a, a) are ascribed 2Na x 2Na gauge matrices, A^°]. We drop in the following the suffix labels (a6), (aa) 
and A, i,j on the matrices, except when needed. We use conventions in which the matrices satisfy the normalization 
and closure sum conditions, Trace(AAA^) = Sab, J2A'^''''^'^^(^^'^A}Trace{02\A} = Trace{Oi02)- 

The orientifold symmetry QTZ acts on the gauge quantum numbers of brane stacks through the twist matrix given 
by the direct product of 2Nfj^ x 2iV^ unitary matrices, ^011,(1, [m = The orientifold projection on physical states is 
then defined by, A^'''' = riA'ynn,a^^A^^'^ Inn b' where rjA are state dependent complex phase factors. These conditions 
must be imposed only for brane stacks n which wrap cycles coinciding with their orientifold mirror images, n = fi' , 
hence fixed under the orientifold twist, 11,^ = 11,^ = Hoe- In the case of a brane stack /i at generic angles, fi 7^ /i', no 
conditions need be imposed beyond requiring that the representations for /i and /i' be conjugate. Nevertheless, it is 
convenient to treat in a unified way the cases with /i = /x' and /x 7^ /i', by taking the orientifold projection matrix in 
the latter subspaces to be trivial, ^Q,n,ii = 1- The SO and Sp orientifold projections, corresponding to Qoe = T2^, 
are characterized by the property of the twist embedding matrix, ^f^^ ^ = ±-fn'R,n- 

We now specialize to the Sp type projection which is the appropriate one for the model of interest to us to be 
discussed in Subsection II B. Solving the condition for the bifundamental representation, A'^'^''^ = —lo.n^n^'^'^'^'''^ Ian w 
with antisymmetric twist matrices, 'ynn,a, [a = of dimension (2A^^+2A^,^) x {2Ni^+2Ni,), and the similar condition 
for the adjoint representation matrix, A^'^''^ of dimension 2A^^ x 2A^^, one obtains 

where the sub-blocks m,a, j3, 7, 5 and si, S2 designate arbitrary generic and symmetric matrices with m, si, S2 
and a, /3, 7, 6 having dimensions A^^^ x N/^ and N/j, x N^, respectively. The conjugate group representations are 
assigned Hermitian conjugate matrices. The adjoint representation matrix A'^'^^'-' in Eq. (A. 10) involves 2Ap(2A^ + l)/2 
independent parameters, as needed to match the dimension of the gauge group USp{2Nfj). The special case for the 
matrix A^''") with si = S2 = ^ corresponds to the adjoint representation of U{N^), involving the expected number of 
A^ independent parameters. For the bifimdamental representation matrix \'^^^\ the entries in the sub-blocks B and 
B' are related by the requirement that the substitution jj ^ jj' corresponds to charge conjugation. The infinitesimal 
transformations of the gauge matrices A^^' in the representation R of sector obtained from the commutator 

with the adjoint representation matrix, 

5^.A(«) = {S^ + 5.)A(«) = [A^'^") e A(«-), A(«)], [A(«-) = (e^, -e^) ® + 7^ (e., -e.)] (A.ll) 

act on the sub-block matrix entries of the matrices B, B' for the bifundamental representations as 

6f_,ya = e^a - ae^, 5^1,(3 = -I- /3e^, (5p^7 = -6^7 - 76,^, = -e^5 -I- St^. (A. 12) 

The four inequivalent bifundamental representations are thus in one-to-one correspondence with the sub-blocks, 
a, (3, 7, 6, characterized by the charge assignments, (Q^, Q,y) : a (1, — 1), /? (1, 1), 7 (—1,-1), 6 (—1, 1). 



APPENDIX B: TWO-BODY PROCESSES AT HIGH ENERGY COLLIDERS 

1. Tree level string amplitudes for processes with fermion and gauge boson pairs 

We discuss in this appendix the tree level four point open string amplitudes in models related to the intersecting 
brane models. For comparison with the results in Subsec. II A, wc first consider the branes within branes models. 
We start, for completeness, with the case of four fermion modes belonging to the diagonal sector {p,p) of a Dp-brane 
which corresponds to the single Dp-hva,ne model [25]. The result can be derived by dimensional reduction of the 
familiar formula for the Z)9-branes of type I theory 

>l(p,p)4 = Gr>pTi234 / dxx~^~'^{l - a;)~*~^[(l - a;)5i234 - a;<Si423] + perms 
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= -Gr,pTi234 J (^«Si234 - s«5i432) + perms, (B.l) 

where we use same notations as in Eqs. (II. 1) and (II. 6). The above formula is formally related to that in Eq. (II. 1) by 
the substitution, C»Si234Vi234 — > —GDp{tSi234 — s<5i432)^^^. Matching to the pole term from massless gauge boson 
exchange determines the normalization factor as, G£)pTi234 = "^Qop' where the £>p-brane gauge coupling constant gjjp 

enters the gauge current vertex as, \/2gDpT°'^\- 

We now focus on the massless fermion modes of the Dp/D{p + 4)-brane models localized at brane intersections in 
the non-diagonal sectors, {p,p + 4) + (p + 4,p). The calculations, discussed initially in [33], bear formal similarities 
with those for intersecting brancs. Wc specialize to the case p = 3 of setups with £'3/D7/-brane pairs with / = 1, 2, 3. 
There are two possible channels for the couplings of four non-diagonal sector modes (3, 7/), which correspond to the 
configurations |(3,7)/|^ and |(3, 7/)p|(3, 7j)p. The string amplitude for the first channel, [(3,7)/]^, involving two 
identical pairs of conjugate fermion modes is given by 



-4(3,7,)^ = C'dp f dxX-'-\l - x)-'-\l - X)S,234 - O^W] (^1234 Ui^^^^J^) + ^^4321 l[{ 

JO \ . r ^ [x) . 



Fi{l-x) 



[ta{x)=z\La\' 'p^^^ M^)=Wj^^y F{x)=Fi-,-;l;x)] (B.2) 

where the label A = {Jm, K,,,). [m = 1,2] in the products runs over the real dimensions of the 4-d sub-torus x T|- 
of the internal torus wrapped by the -D7/-brane, the Wilson line parameters along the corresponding sub-torus are 
denoted by ca, and the sub-torus volume parameter is defined by, ILa]"^ = r'jr\. The direct and reverse orientation 
terms inside the large parentheses arc related by the change of integration variable, x (1 — .t). Combining the regions 
of the x-integral near x = Q and x = 1 yields the low-energy approximate representation of the string amplitude as 
infinite series of s-channel or f-channel poles located at the string compactification modes 

T'l234[»5l234 , , ; To-o «54321 . , Tl~2 ] 

p^e^V -s + Y.AiPA + eA?r\ -t + Y.APAlrA 

+I432l['5i234 — 2 / 2 '^4321 , , v-^ / j \o 2 ] 1 ' (^-3) 

where the momentum modes in the open string sector {7i,7i) arise after use of the Poisson resummation formula 
and the winding modes belong to the open string sector (3,3). The massless pole terms determine the normalization 

constant in terms of the gauge coupling constants by the formula, C^pTi234 = ^ngs = ^^°(2ot)p-™° — with Vp-3 ~ 
(27rr)^'~^ denoting the volume of the {p — 3)-cycle of the internal manifold wrapped by the Dp-brane. For the Abelian 
gauge group case, C^pTi234 = C^p?432i = 2g%^ = 2g%^r'^jr%. 

For the channel |(3, 7/)p|(3, 7j)p, involving two distinct pairs of conjugate fermion modes, (3, 7/) and (3, 7j), [I ^ 
J] the string amplitude is given by the formula 



■^(3,7,)^(3,7.)= = C'op I dxX-'-\l - (ti234[(1 - x)5i234 + xVu32\ 

!(1-X) 



(^) 

+T432i[(l " a;)Pi234 + x5i432] H^M^)) ' (^-4) 

F2(i-x) 1 



where 



^^1234 = («r7°«2)(«37°«4), 5i234 = («?^7°7''^^2)(«^7%m«4). (B.5) 

The label A in the products Ha '^^'^^ o^^^^ directions of the wrapped internal sub-torus common to the 
i)7/, _D7.7-branes, corresponding to ^4 = \<n = 1,2] and ta, tb retain the same definition as above except 

that the volume parameter is now given by, = t'k- The Dirac spinor scalar quartic coupling, 7^1234, appears 

because the modes in the non-diagonal open string sector (7/,7j) to which the fermion pairs couple are Lorentz 
scalars. Combining the contributions from the regions of the a;-integral near x = and x = 1, yields the low-energy 
representations as infinite series of s-channel and f-channel poles located at the string compactification modes 

-4'(3,7.)^0 + -4'(3,7.)M - C'j,^ 1^ (^Ti234[5r234_^^^^(^^^^^^2,2 + ^"^^ -t + Y.AP\lr\ ^ 
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Y.AV\lr\ ^ -t + EAiPA + CAn 



+74321 ['P1234 , — ir7l2 '^1432 + , f„ , I , \2 2 ] ) • 

(B.6) 



The normalization factor is given by the same formula as found above. 

We next discuss the 2^2 body processes involving the gauge boson pair production by fermion-antifermion 
annihilation and the gauge boson pair scattering processes. The string amplitudes have a universal form with the 
model dependence residing only in the gauge structure. In particular, identical formulas hold in single and multiple 
Dp-hT&ne models. For the localized fermions of intersecting brane models, the string amplitude for photon pair 
production by fermion-antifermion annihilation, A'^^^^ = ^'(e+(fci) + e~{k2) + 7(^3) + 7(^4)), is calculated from 

the world sheet vacuum correlator, < V^g ^''Vg ^'v^^ ^'V"!"' >. The dependence on the interbranc angles from the 
correlator of a single pair of coordinate twist fields is found to exactly cancel that coming from contracting the 
corresponding pair of spinor twist fields. The rest of the calculation is standard and gives the same result as that 
obtained by dimensional reduction from the £>9-brane amplitude [91] 

A' -r IT '^(^'^) T ^CM) ^(fi^i^r^, . . \ 

■^{ee-y-i) " ^Dp[-'l234 -'1324 7 T -il243 J-f^ (,"1 , "2, Cs, e4j, 

Sb U/O sill 

[K{ui,U2, 63, 64) = t {Ui^sih + ^i)h'^2) + U (Uie4(^2 + ^3)j^3W2), 

T1234 = Trace(AiA2A3A4) + Troce(A4A3A2Ai)] (B.7) 

where the factor depending on the polarization wave functions, K{ui , U2 , £3 , £4) , is (anti) symmetric imder permutations 
of the (fermion) boson particle labels. Assuming, for simplicity, the three gauge trace factors to be equal, the 
factorization on the massless pole terms identifies the normalization constant to the gauge coupling constant in the 
Abehan and U{N) non-Abelian group cases as, GDpTi2iA = 2g%p and GDpTi234 = 25|5p[Ea(^")i2- (^")34 + ^(l)i2- 
(1)34]. 

The string amplitude for the gauge boson pair scattering process, A'^^^^ = A'{'-f{ki) +7(^2) +7(^3) +7(^4)), is of 
same form as that obtained in the familiar £>9-brane case [91] 

^(^^^^) = G'jjp[Ti234-^S{s,t) +Tis24-^S{u,t) + Ti243—S{u, s)]K^{e-i,e2,e3,e4), 

SL UjL SU 

1 s 
K^{ei,e2, 63, €4,) = --^[st (ei • e3)(e2 • €4) + perms] + - ((d • fc4)(e3 • A:2)(e2 • €4) + perms), (B.8) 

where the normalization constant is related to the gauge coupling constant by, G^p7i234 = G^p7i324 = G'jjpTi243 = 

^ghp- 



2. Helicity amplitudes 

To enable the comparison with experimental measurements, it is useful to express the various string ampli- 
tudes in the spin helicity basis of the various modes. A convenient way to proceed is by first establishing the 
correspondence dictionary between the kinematics of string amplitudes, where all particles are incoming, with 
that of physical processes, and using next the familiar crossing relations which transform particles to antiparti- 
cles and flip the sign of momenta and helicities. Thus, the amplitude for fermion-antifermion pair production, 
fiipi) + f2iP2) /.-rfe) + ft{Pi)^ is obtained from the string amplitude, A{Ji{ki) + /2(fc2) + f?,{k3) + fi{k4)), by 
setting, /i(fci) = ft{pi), f2{k2) = /2"(P2), 73(^3) = /"(-P3), /4(A;4) = /+(-P4), which involves substituting the 
momenta, kinematic variables and Dirac spinors as 



[fcl, k2, fcs, k4] =^ [pi, P2, -P3, -P4], 

[s = -(fcl + fcaf , t = -(fc2 + ksf, u = 

[Ul(fcl), U2(fc2), ■"3(^3), W4(fc4)] => [V*{pi), U{P2), U*{pz), v{p4)\- (B.9) 



|S = -(fci + fc2f, ^ = -(fc2 + fc3)^ M=-(fcl + fc3)^] =^ [S = -(P1+P2)^ ^ = -(Pl-p4)^ U = -{p-l - psf], 



For the choice of kinematic variables in the center of mass frame, pi = —p2 = p, —p?, = P4 = jf , [p ■ p' = cos6] the 
kinematic invariants read, s = ^p^, t = — ssin^ |, u = — scos^ |. The helicity polarization basis for the spin one- half 
Dirac fermions is described by the familiar formulas 

e = e" : u{p, A) = -pQ^ ® <t>x{p); e = e+ : v{p, A) = ig) ^-a(p), 

(— sin-\ /cos-\ 
cos I / Vsin|y 
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The Dirac spinor matrix elements in the configuration of helicities for the physical process, /+/ —>//"*" are given 

by 





(+-- + ) 


(+-+-) 


( ) 




-2u 


-2t 





<Si432 = {v{pi)'y^v{p4)){u{p3)'yi^u{p2)) 


2u 





-2s 


i«5l234 — S«5l432 


2^2 


-2*2 


2s2 



The results for the other helicity configurations are inferred from the above formulas by invoking the symmetry 
under space parity. The string amplitude in the three independent helicity configurations for the physical process, 
cai + ~* + eAg , are given in the single brane and the intersecting brane cases by 



mTp> = -Gdp^^[2u', -2t\2s'] + perms, 

(^'/4+" + "'+""+' ')/SS = 2^5,Ti234Vi234[-2u, -2t, 0]. (B.ll) 



To account for the electroweak symmetry breaking in the charge neutral channels one needs to substitute the massless 
photon pole term by the sum of 7 + Z boson pole terms. The Z boson exchange contribution is obtained from that 
of 7 exchange by the substitution, s ^ s — m%, along with the following replacements for the chirality couplings 
LL^ RR^ LR^ RL'. 

e} ^ e}[aL{.f)alif), ai,if)a*j,{f), aL{f)a*M), afl(/)al(/)], 

_i _l_ ^,2 

[aL{f) = — —, = — , sw = sm9w,cw = cos9w]- (B.12) 

SwCw Cw 

The diflcrcntial cross section for the spin-impolarized Bhabha scattering process. e + e e + e, obtained by adding 
to the Standard Model terms the contributions from the 2? = 6 contact interactions, is given by 

All = A{e+el ^ e+el) = A+^+- + = ] + + 

Arr = A{e+e^ -> e+e^) = As+-+ + [As+-+ = ^ + j + ' ^ 



t ■ 



Arl = A{e+e^ ^ e+e^) = Asm +^,[Asm = At^++ = ^ + aL{e)aR{e) ^ 



Arl = A{e^ej^ ^l^r) = ^sm + ' ^"^^m = ^sm^ = ~ + «ii(e)aL(e)^ - m?'^' 

At - A(r+r- ; r+r^] - ,4++++ I 
^LR — ^K'^r'^^r ^ '^^r'^-rJ — ^SM + ^^2 ' 

Alr = A{e\ej^—> e^e^) = Ag^ oA^' (B.13) 

where a = I7 and denote the coefficients of the local operators previously defined in Eq. (11.17), and anie) 

designate the Z-boson vertex couplings defined in Eq. (B.12). We use the suffix label s, t to distinguish the s- and 
i-channel pole terms. 

For the fermion pair scattering processes, the correspondence dictionary between the kinematical variables in the 
physical process, A{fi{pi) + h{P2) hiPs) + Mpa)), and the string process, A{fi{ki) + 72(^2) + 73(^3) + fi{ki)), 
can be written as 

[ki,k2,ks,ki]^[-ps,Pi,-p4„P'2], [wi(A;i), U2(A;2), -"3(^:3), ^4(^4)] =^ [■«*(P3), u{pi), u*{pa), m(P2)], 

[s = -(fcl+fc2)^ i=-{k2 + ksf, U=-{ki+ksf\^[t = -{j>i-pif, u=-{pi-p3f, s = -{pi + 1(0^1.4) 

For clarity, we have distinguished the kinematic invariants of the string theory process by adding momentarily hat 
symbols. In the center of mass frame of the physical process, pi = —p2 = P, —Pz — Pi — p' , [p ' p' = cosO] 
the kinematic invariants read, s — Ap^, t — — ssin^ |, u = — gcos^ |. We specialize now to the physical process, 
&h{pi)+<1h' {P2) — > eH{p'i)+qH' [Pi), where the Dirac spinor matrix elements in the single D-brane and the intersecting 
brane cases, signalled by the suffix labels Dp and ISB, are given by 

(i5i234 - s5i432)dp = ps^, -2u^]; (5i234)/SB = [-2u, -2i] = [-2s, -2u], (B.15) 
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with the two entries inside brackets corresponding to the equal and unequal helicity cases, [H = H' , H ^ H']. The 
helicity basis amplitudes for the physical process can be written for the single D-brane and the intersecting brane 
cases as 

s s f 1 \ 

■A{eH + QH' ^ Sh + qH')Dp = GDp[S{u,t)—Ti324+S{s,u)-Ti234 + S{t,s)-Ti243][ _ 2 / 2 )' 

A{eH + qw BH + qH')lSB = -C[Vl324(w,t)Ti324 + Vl234(s,u)Ti234 + Vl243(s,t)'?"l243]( ) , (B.16) 



2uJ' 

with the upper and lower entries corresponding to the configurations with equal and unequal helicities, [H ^ H' , H ^ 
H']. 

The helicity amplitudes for the gauge boson pair production process, e+(pi) + e (^2) 7(^3) + 7(^4)) 8i,ie given 
by the formulas 

-4(4 + eZ ^ 7i +7ii) = -jAeR + el ^ 1r + 1l) = 2gl^^ -[-S{s,t) - S{u,t) + -S{s,u)], (B.17) 

where we have assumed the various gauge factors to be equal. The helicity amplitudes for the gauge boson pair 
scattering process are given by the explicit formulas 

[/(s, t, u) = sS{t, u) + tS{u, s) + uS{s, t)] (B.18) 

where the upper and lower entries correspond to the configurations with equal and unequal helicities, [H = H' , H ^ 
H'\. Note the crossing relations, ^(7il + 72/? 73L + HR)\s,t,u = A{-iiL + J4R 73L + l2R)\t,s,u- 



[1] E.G. Gimon and J. Polchinski, Phys. Rev. D54, 1667 (1996) [arXiv:hep-th/9601038] 

[2] E.G. Gimon and C.V. Johnson, Nucl. Phys. B477, 715 (1996) [arXiv:hep-th/9604129]; Nucl. Phys. B479, 285 (1996) 
[3] M.R. Douglas and G. Moore, [arXiv:hep-th/9603167] 

[4] M. Berkooz, M.R. Douglas, and R.G. Leigh, Nucl. Phys. B480, 265 (1996); V. Balasubramanian and R.G. Leigh, Phys. Rev. 
D55, 6415 (1997) 

[5] E. Dudas, Class. Quant. Grav. 172000R41-R116 [arXiv:hep-ph/0006190]; C. Angelantonj and A. Sagnotti, Phys. Rep. 

371, 1 (2002), Erratum-ibid. Phys. Rep. 376, 339 (2003) [arXiv:hcp-th/0204089] 
[6] R. Blumenhagen, M. Cvetic, P. Langackcr, and G. Shiu, [arXiv:hep-th/0502005] 
[7] R. Blumenhagen, B. Kors, D. Lust, and S. Stieberger, [arXiv:hep-th/0610327] 
[8] F. Marchesano, [arXiv:hep-th/0702094] 
[9] I. Antoniadis, Phys. Lett. B246, 377 (1990) 
[10] J. Lykkcn, Phys. Rev. D54, 3693, (1996) [arXiv:hep-ph/9603133]; E. Witten, Nucl. Phys. B471, 135 (1996) [arXiv:hep- 

th/9602070] 

[11] N. Arkani-Hamcd, S. Dimopoulos, and G. Dvah, Phys. Rev. D59, 086004 (1999) [arXiv:hep-ph/9807344]; I. Antoniadis, 

N. Arkani-Hamed, S. Dimopoulos, and G. DvaU, Phys. Lett. B436, 257 (1998) [arXiv:hep-ph/9804398] 
[12] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999); Phys. Rev. Lett. 83, 4690 (1999) 
[13] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D61, 033005 (2000) [arXiv:hep-th/9903417] 
[14] N. Arkani-Hamed, Y. Grossman, and M. Schmahz, Phys. Rev. D61, 115004 (2000) 
[15] E.A. MirabeUi and M. Schmaltz, Phys. Rev. D61, 113011 (2000) [arXiv:hep-th/9912265] 
[16] H. Georgi, A.K. Grant, and G. Hailu, Phys. Rev. D63, 064027 (2001) 

[17] T. Han, J.D. Lykkcn, and R.J. Zhang, Phys. Rev. D59, 105006 (1999) [arXiv:hep-ph/9811350] 
[18] G.F. Giudice, R. Rattazzi, and J.D. Wells, Nucl. Phys. B544, 3 (1999) [arXiv:hep-ph/9811291] 
[19] J.L. Hewett and M. Spiropulu, Ann. Rev. Nucl. Part. Sci. 52, 397 (2002) [arXiv:hep-ph/0205106] 

[20] K. Cheung, [arXiv:hep-ph/0305003] 

[21] G. Landsberg, SLAG Summer Institute, 2004, [arXiv:hep-ex/0412028] 

[22] Z. Kakushadzo, Nucl. Phys. B552, 3 (1999) [arXiv:hop-th/9812163]; Nucl. Phys. B548, 205 (1999) [arXiv;hep-th/9811193]; 

Z. Kakushadze and S.H. Henry Tye, Nucl. Phys. B548, 180 (1999) [arXiv:hep-th/9809147]; Z. Kakushadze and S.H. Henry 

Tye, Phys. Rev. D58, 126001 (1998) [arXiv:hep-th/9806143] 
[23] E. Kiritsis, Fortsch. Phys. 52, 200 (2004); Phys. Rep. 421, 105 (2005) [arXiv:hep-th/0310001] 
[24] E. Dudas and J. Mourad, Nucl. Phys. B575, 3 (2000) [arXiv:hep-th/9911019] 

[25] S. Cullon, M. Porolstein and M.E. Peskin, Phys. Rev. D62, 055012 (2000) [arXiv:hop-ph/0001166] 
[26] F. Cornet, J.I. Illana, and M. Masip, Phys. Rev. Lett. 86, 4235 (2001) [arXiv:hep-ph/0102065] 



27 



J. Frioss, T. Han and D. Hooper, [arXiv:hep-ph/0204112] 

P. Burikham, T. Han, F. Hussain, and D.W. McKay, [arXiv:hcp-ph/0309132] 

T. Burikham, T. Figy, and T. Han, [arXiv:hcp-ph/0411094] 

P. Burikham, [arXiv:hep-ph/0502102]; [arXiv:hep-ph/0601142] 

L.A Anchordoqui, H. Goldberg, D. Lust, S. Nawata, S. Stieberger, and T. R. Taylor, [arXiv:hep-ph/0808.0497] 

D. Lust, S. Stieberger, and T.R. Taylor, [arXiv:hep-ph/0807.3333] 

L Antoniadis, K. Benakh, and A. Laugier, JHEP C0105, 044 (2001) [arXiv:hcp-th/0011281] 
M. Bershadsky and A. Radul, Int. J. Mod. Phys. A2, 165 (1987) 

L.J. Dixon, D. Fricdan, E.J. Martinec, and S.H. Shcnkcr, Nucl. Phys. B282, 13 (1987) 
A. Hashimoto, Nucl. Phys. B496, 243 (1997) [arXiv:hep-th/9608127] 

E. Cava, K.S. Narain, and M.H. Sarmadi, Nucl. Phys. B504, 214 (1997) [arXiv:hep-th/9704006l 

J. Frolich, O. Grandjean, A. Recknagel, and V. Schomerus, Nucl. Phys. B583, 381 (2000) [arXiv:hep-th/9912079] 

M. Cvetic and L Papadimitriou, Phys. Rev. D68, 046001 (2003) [arXiv:hep-th/0303083]; M. Cvetic and I. Papadimitriou, 

Phys. Rev. D79, 029903 (2004) (Erratum) 

S.A. Abel and A.W. Owen, [arXiv:liep-th/0205031] 

S.A. Abel and A.W. Owen, Nucl. Phys. B663, 197 (2003) [arXiv:hep-th/0303124] 
S.A. Abel and A.W. Owen, [arXiv:hep-th/0310257] 

I.R. Klebanov and E. Witten, Nucl. Phys. B664, 3 (2003) [arXiv:hep-th/0304079l 
M. Chemtob, [arXiv:hep-ph/0702065] 

F. Cachazo and P. Svrcck, [arXiv:hep-th/0504194] 

S. Stieberger and T.R. Taylor, [arXiv:hep-th/0708.0574] 

I. Antoniadis, [arXiv:hep-th/0102202]; I. Antoniadis, K. Benakli, Int. J. Mod. Phys. AA15, 4237 (2000) [arXiv:hep- 
ph/0007226] 

G. Shiu, R. Schrock, and S.H.H. Tye, [arXiv:hep-ph/9904262] 

S. Abel, M. Masip, J. Santiago, JHEP C0304, 057 (2003) [arXiv:hep-ph/0303087] 

S. Abel, O. Lebedev, and J. Santiago, Nucl. Phys. B696, 141 (2004) [arXiv:hep-ph/0312157] 

A. Delgado, A. Pomarol, and M. Quiros, [arXiv:hep-ph/9911252] 

T.G. Rizzo, Phys. Rev. D64, 015003 (2001) [arXiv:hep-ph/0101278] 

Y. Grossman and G. Perez, Phys. Rev. D67, 015011 (2003) [arXiv:hcp-th/0210053] 

A. Masiero, M. Peloso, L. Sorbo, and R. Tabbash, Phys. Rev. D62, 063515 (2000); D.E. Kaplan and T.M. Tait, JHEP Cll, 
051 (2001) [arXiv:hep-ph/0110126]; S. Nussinov and R. Schrock, Phys. Lett. B526, 137 (2002); S. Nussinov and R. Schrock, 

Phys. Rev. Lett. 88, 171601 (2002); D.J.H. Chung and T. Dent, Phys. Rev. D66, 023501 (2002) [arXiv:hep-ph/0112360] 
G. Barcnboim, G.C. Branco, A. dc Gouvca, and M.N. Rebelo, Phys. Rev. D64, 073005 (2001) [arXiv:hep-th/0104312]; 
G.C. Branco, A. dc Gouvca, and M.N. Rcbclo, Phys. Lett. B506, 115 (2001) 

B. Lillie and J.L. Hewett, Phys. Rev. D68, 116002 (2003) [arXiv:hep-ph/0306193] 

B. Lillie, JHEP C12, 030 (2003) [arXiv:hep-ph/0308091] 

G. Moreau and J.L Silva-Marcos, JHEP C0603, 090 (2006) [arXiv:hep-ph/0602155] 

C. Csaki, A. Falkowski, and A. Weiler, JHEP C09, 008 (2008) [arXiv:hep-ph/0804.1954]; [arXiv:hep-ph/0806.3757] 

D. Crcmadcs, L.E. Ibancz, and F. Marchcsano, [arXiv:hep-ph/0212064] 

D. Cremades, L.E. Ibanez, and F. Marchcsano, JHEP C0307, 038 (2003) [arXiv:hep-th/0302105] 
C. Kokorelis, [arXiv:hep-th/0309070] 

M. Cvetic, T. Li and T. Liu, Nucl. Phys. B698, 163 (2004) [arXiv:hep-th/0403061] 

F. Gmeiner, R. Blumenhagen, G. Honecker, D. Lust, and T. Weigand, JHEP COl, 004 (2006) [arXiv:hep-th/0510170] 
M.R. Douglas and W. Taylor, JHEP COl, 031 (2007) [arXiv:hep-th/0606109] 

F. Grncincr and G. Honecker, [arXiv:hcp-th/0806.3039] 

T. Higaki, N. Kitazawa, T. Kobayashi, and K.-J. Takahashi, Phys. Rev. D72, 086003 (2005) [arXiv:hep-th/0504019] 

T.T. Burwick, R.K. Kaiser, and H.F. Miiller, Nucl. Phys. B355, 689 (1991) 

J. Polchinski, 'String Theory' Vols. I, II (Cambridge University Press, Cambridge, 1998) 

M. Bando, T. Kugo, T. Noguchi, and K. Yoshioka, Phys. Rev. Lett. 83, 3601 (1999) [arXiv:hep-ph/9906549] 

J. Hisano and N. Okada, Phys. Rev. D61, 106003 (2000) [arXiv:hep-ph/9909555] 

H. Murayama and J.D. Wells, Phys. Rev. D65, 056011 (2002) [arXiv:hep-ph/0109004] 

E. J. Eichten, K.D. Lane, and M.E. Peskin, Phys. Rev. Lett. 50, 811 (1983) 

N. Chamoun, S. Khalil, and E. Lashin, Phys. Rev. D69, 095011 (2004) [arXiv:hep-ph/0309169] 

R. Blumenhagen, D. Liist, and S. Stieberger, [arXiv:hep-th/0305146]; D. Liist, [arXiv:hep-th/0401156] 

P. Anastasopoulos, T.P.T. Dijsktra, E. Kiritsis, and A.N. Schellekens, Nucl. Phys. B759, 83 (2006) [arXiv:hep-th/0605226] 

K.R. Dienes, E. Dudas, and T. Gherghetta, Phys. Lett. B436, 55 (1998) [arXiv:hep-ph/9803466]; Nucl. Phys. B537, 47 

(1999) [arXiv:hep-ph/9806292] 

Joanne L. Hewett, Phys. Rev. Lett. 82, 4765 (1999) 

T.G. Rizzo, Phys. Rev. D59, 115010 (1999) [arXiv:hep-ph/9901209] 

I. Antoniadis, K. Benakh, and M. Quiros, Phys. Lett. B460, 176 (1999) [arXiv:hep-ph/9905311] 
Particle Data Group Collaboration, Jour. Part. Phys. 33, 1155 (2006) 

J.A. Grifols, E. Masso, and R. Toldra, Phys. Rev. D57, 2005 (1998) [arXiv:hep-ph/9707531] 
K. Cheung and Y.-F. Liu, Phys. Rev. D72, 015010 (2005) [arXiv:hcp-ph/0505241] 

The LEP CoUaborations: ALEPH Collaboration, DELPHI Collaboration, L3 Collaboration, OPAL CoUaboration, the 



28 




FIG. 1: The contact interaction coefficient g^^i = mlKij_/{PJ^^,g'^) for the operator e\ predicted in Cremades et ai, [60, 61] 
model is plotted as a function of mgr. Same results hold for the operator q% and numerically close results hold for the operators 
u%, d% and e%. The reference set of parameters specifying the model is defined by, p = 1, e = l,e = 1, ~ li Pi = 1) P2 = 
1, CA = (0,0,0), es — (0,0,0) (see the text just below Eq. (ILSO)). In panel (a), the five curves I , II , III , IV,V are obtained 
by starting from the reference set (I) and performing in succession the following variations: = 2 {II), P2 ~ 1/2 (III), 
/3i = /32 = 1/2 {IV) and x' ~ ^ {^)- ^'^ panel {b), the three curves I, II, III are obtained by starting from the reference set of 
parameters (/) and including the changes, es = (0, i,0) (77) and en = (0, |,0) {III). In panel (c), the two curves I, II refer 
to the reference set of parameters with the choices, eA = ~ (0, |,0) (7) and eA ~ es = (0, |,0) {II). The presence of cusps 
in the curves of panel (c) is due to the use of semi-logarithmic plots for the absolute values of the coefficients. 
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FIG. 2: The contact interaction coefficient gfj^i for the operator (e^)* ~ e\ is plotted as a function ofrusr. We use same 
conventions as in Fig. (1). Same results hold for the operator (d^)^ ~ d%^ and numerically close results hold for the operator 
(u^)* ~ The three panels (a), (&), (c) are associated to = |j Ij 2. In each panel, we display three curves obtained by 
starting from the reference set of parameters (/) and changing the longitudinal distance parameter to cb = (0, ^,0) {II) and 
cb ~ (0, |, 0) {III). The presence of cusps in certain curves is due to the use of semi-logarithmic plots for the absolute values 
of the coefficients. 
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(a) =e3=(0,0,0) (b) e^^'e^', £^=(0,0,0) (c) f^^j^X'' =h=(^fi'^) 




FIG. 3: The contact interaction coefficients gfj^i ~ m^gKfj_/{K^J^,g'^) are plotted as a function ofnisr for various configurations 
of the quarks and leptons. We use same conventions as in Fig. (1). In the panel (a) referring to the operator e|,(e^)^ ~ eie/j, 
we display four curves obtained by starting from the reference set of parameters (I), and independently changing the single 
parameters, = 2 (//), 02 = 1/2 (HI), (3i = (32 = 1/2 {IV). In the panel (6) referring to the operator e\e\, we display 
two curves obtained for the reference set of parameters with the values of the longitudinal distance parameter, es = (0, ^,0) (/) 
and es = (0, §,0) (//)• In panel (c), we display the coefficients for the quark and quark-lepton operators, q^u\, u\e%, d^e^, 
obtained with the reference set of parameters. 
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FIG. 4: The ratio (da / d cos 9) s M+contact/ {da /d cos 6) sm for the Bhabha scattering differential cross section at the center of 
mass energy ^/s — 183 GeV is plotted as a function of cos 6 using the coefficients of contact interactions predicted in Cremades 
et al, [60, 61] with the reference set of parameters. We consider four values of the string scale parameter, rus = 0.5, 1.5, 2., 3. 
TeV, and two values of the compactification scale parameter, rugr = 1, 3. The group of four lowermost curves from the bottom 
right comer is associated to rusr = 3, and the group of four uppermost curves close to the horizontal axis is associated to 
rusr = 1. The predictions for variable string scale within each group are drawn with dashed curves using dashings of increasing 
length in correspondence with the increasing sequence of values, rUs = 0.5, 1.5, 2., 3. TeV. 
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FIG. 5: The ratio (da / d cos 9) s M+contact/ {da /d cos 0)sm for Bhabha scattering differential cross section at the center of mass 
energy ^/s = 500 Ge V is plotted as a function of cos 6 using the coeffi,cients of contact interactions predicted in Cremades et 
al, [60, 61] with the reference set of parameters. The group of four lowermost curves from the bottom right comer is associated 
to msV = 3, and the group of four uppermost curves close to the horizontal axis is associated to msV = 1. The predictions for 
variable string scale within each group are drawn with dashed curves using dashings of increasing length in correspondence with 
the increasing sequence of values, rus = 2., 3., 4., 5. TeV. 
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FIG. 6: The ratio {da /d cos 0)sM+contact/ {da /d cos 9) sm for Bhabha scattering differential cross section at the center of mass 
energies ^/s — 183. GeV and = 500. GeV are plotted as a function of cos 8 for the string amplitudes predicted in Cremades 
et ai, [60, 61] model with the reference set of parameters, using g^ = g2{'mz) = 0.425 and a{mz) ~ 1/127.9. The total cross 
section is evaluated for the subtraction regularized string amplitudes with the massless gauge boson pole terms replaced by the 
corresponding pole terms at the physical masses of the neutral 7, Z gauge bosons. We consider the approximate estimate for 
the mixed chirality amplitudes, G^^, G^^ , setting these to a constant multiple of the pure chirality amplitudes, G^^ , G^^ , 
defined by the parameterization, G^^ = G^^ = xG^^ = xG^^, with the two extreme numerical values x = ^, 5. We 
also consider the two values for the string mass scale, tUs — 1, 2 TeV and rris = 2, 4 TeV, in correspondence with the two 
center of mass energies, y/s = 0.183 TeV and y/s = 0.500 TeV. On the left hand side, the upper and lower panels (a) and 
(6) display the ratio for mar = 1 and rrisr — 3 at ^/s = 183. GeV with the four curves I, II, III, IV (in full, dotted, 
short-dashed, dashed lines) referring to the values of the string mass scale and the proportionality factor between the mixed and 
pure chirality amplitudes: {x = G^^ /G^^ ,ms/ TeV) = (|, 1), (^,2), (5, 1), (5,2) TeV. On the right hand side, the upper and 
lower panels (c) and (d) display the ratio for nisr — 1 and m^r = 3 at ^/s — 500. GeV with the four curves I, II, III, IV 
referring to the values of the string mass scale and the proportionality factor between the mixed and pure chirality amplitudes: 
(a; = G^«/G",m,/reV) = (i,2), (1,4), (5,2), (5,4). 
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FIG. 7: Lower bounds on rus as a function of maV deduced from experimental data for the neutral mesons mass shifts and the 
three-body leptomc decays of charged leptons. The curves labeled I, II, III, IV refer to the \ AFq\ — 2 mass shift observables 
of the neutral meson systems K — K , B — B , Ba — Bs, D — D using the experimental inputs quoted in Eq, (IV. 11), The curves 
labeled V, VI refer to the |A_F;| = 1 observables for the charged leptons three-body decay rates, /i^e + e + e, r^e + e + e 
using the experimental inputs quoted in Eq. (IV. 12). 



